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Preface 



Gravitational waves and their detection are becoming more and more im- 
portant both for the theoretical physicist and the astrophysicist. In fact, 
technological developments have enabled the construction such sensitive 
detectors (bars and interferometers) that the detection of gravitational ra- 
diation could become a reality during the next few years. In these lectures 
we give a brief overview of this interesting and challenging field of modern 
physics. 

The topics to be covered are divided into 6 Lectures. We begin by de- 
scribing gravitational waves in linearized general relativity, where one can 
examine most of the basic properties of gravitational radiation itself: prop- 
agation, gauge invariance, and interactions with matter (and in particular 
with detectors). 

The second Lecture deals with gravitational waves detectors: how they 
operate, what their most important sources of noise are, and what mech- 
anisms are used to overcome noise. We report here on the most impor- 
tant detectors planned or under construction (both ground-based ones and 
space-based), their likely sensitivity, and their prospects for making detec- 
tions. Other speakers will go into much more detail on specific detectors, 
such as LISA. 

The third Lecture deals with the astrophysics of likely sources of grav- 
itational waves: binary systems, neutron stars, pulsars. X-ray sources, su- 
pernovae/hypernovae, 7-ray bursts and the big bang. We estimate there 
the expected wave amplitude h and the suitability of specific detectors for 
seeing waves from each source. 

The fourth Lecture is much more theoretical. Here we develop the 
mathematical theory of gravitational waves in general, their effective stress- 
energy tensor, the energy carried by gravitational waves, and the energy in 
a random wave field (gravitational background generated by the big bang). 

The fifth Lecture carries the theory further and examines the gen- 
eration of gravitational radiation in linearized theory. We show in some 
detail how both mass-quadrupole and current-quadrupole radiation is gen- 
erated, including how characteristics of the radiation like its polarization 
are related to the motion of the source. Current-quadrupole radiation has 
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become important very recently and may indeed be one of the first forms 
of gravitational radiation to be detected. We attempt to give a physical 
description of the way it is generated. 

The final Lecture will explore applications of the theory we have de- 
veloped to some sources. We calculate the quadrupole moment of a binary 
system, the energy radiated in the Newtonian approximation and the back- 
reaction on the orbit. We conclude with a brief introduction the current- 
quadrupole-driven instability in the r-modes of neutron stars. 

Chapters one and four are followed by a few exercises to assist students. 
We presume the reader has some background in general relativity and its 
mathematical tools in differential geometry, at the level of the introductory 
chapters of Schutz (1985). See the list of references at the end of these 
lectures for sources suitable for further and background reading. 



B.F. Schutz, F. Ricci 

May 2000 



Chapter 1 

Elements of gravitational waves 



General relativity is a theory of gravity that is consistent with special rel- 
ativity in many respects, and in particular with the principle that nothing 
travels faster than light. This means that changes in the gravitational 
field cannot be felt everywhere instantaneously: they must propagate. In 
general relativity they propagate at exactly the same speed as vacuum 
electromagnetic waves: the speed of light. These propagating changes are 
called gravitational waves. 

However, general relativity is a non-linear theory and there is, in gen- 
eral, no sharp distinction between the part of the nietic that represents 
the waves and rest of the metric. Only in certain approximations can we 
clearly define gravitational radiation. Three interesting approximations in 
which it is possible to make this distinction are: 

• linearized theory; 

• small perturbations of a smooth, time-independent background metric; 

• post-newtonian theory. 

The simplest starting point for our discussion is certainly linearized 
theory, which is a weak-field approximation to general relativity, where the 
equations are written and solved in a nearly flat space-time. The static and 
wave parts of the field cleanly separate. We idealize gravitational waves as 
a "ripple" propagating through a flat and empty universe. 

This picture is a simple case of the more general "short-wave approx- 
imation" , in which waves appear as small perturbations of a smooth back- 
ground that is time-dependent and whose radius of curvature is much larger 
than the wavelength of the waves. We will describe this in detail in Lec- 
ture 4. This approximation describes wave propagation well, but it is in- 
adequate for wave generation. The most useful approximation for sources 
is the post-Newtonian approximation, where waves arise at a high order in 
corrections that carry general relativity away from its Newtonian limit; we 
treat these in Lectures 5 and 6. 

1 



Mathematics of linearized theory 2 

For now wc concentrate our attention on linearized theory. We follow 
the notation and conventions of Misner, et al, (1973) and Schutz (1985). In 
particular we choose units in which c = G = 1; Greek indices run from to 
3; Latin indices run from 1 to 3; repeated indices are summed; commas in 
subscripts or superscripts denote partial derivatives; and semicolons denote 
covariant derivatives. The metric has positive signature. See these two 
textbooks or others referred to at the end of these lectures for more details 
on the theory that we sketch here. For an even simpler introduction, based 
on a scalar analogy to general relativity, see reference [Tj. 

1.1 Mathematics of linearized theory 

Consider a perturbed flat space-time. Its metric tensor can be written as 

gal3 =rial3 +hal3, | /la/3 | < 1 a,/3 = 0,...,3 (1.1) 

where r]ai3 is the Minkowski metric (-1,1,1,1) and hap is a very small pertur- 
bation of the flat space-time metric. Linearized theory is an approximation 
to general relativity that is correct to first order in the size of this perturba- 
tion. Since the size of tensor components depends on coordinates, one must 
be careful with such a definition. What we require for linearized theory to 
be valid is that there should exist a coordinate system in which Equation 
p. II) holds in a suitably large region of space-time. Even though riap is not 
the true metric tensor, we are free to define raising and lowering indices 
of the perturbation with rja/j, as if it wree a tensor on flat space-time. We 
write 

This leads to the following equation for the inverse metric, correct to first 
order (all we want in linearized theory): 

^a/3 ^ ^a/3 _ ^a/3_ ^^3) 

The mathematics is simpler if we define the trace-reversed metric per- 
turbation: 

hap ■■= ha/3 - Tflalih, (1-3) 

where h :— rjaphf^^ ■ There is considerable coordinate freedom in the com- 
ponents hap , since we can wiggle and stretch the coordinate system with a 
comparable amplitude and change the components. This coordinate free- 
dom is called gauge freedom, by analogy with electromagnetism. We use 
this freedom to enforce the Lorentz (or Hilbert) gauge: 

h°'^,l3 = (1.4) 



Using TT gauge to understand gravitational waves 3 

In this gauge the Einstein field equations (neglecting the quadratic and 
higher terms in ft,"^) are just a set of decoupled linear wave equations: 

l^+V^^/i"^ =-16^T"'5. (1.5) 

ot^ J 

To understand wave propagation we look for the easiest solution of the 
vacuum gravitational field equations: 

nh''^=f-^+\7Ah''^^0 (1.6) 

Plane wave solutions have the form: 

hai3 — Asais exp(ifc^x^), (1.7) 

where the amplitude A, polarization tensor e"^ and wave vector k^ are all 
constants. (As usual one has to take the real part of this expression.) 

The Einstein equations imply that the wave vector is "light-like", 
k'^k^ = 0, and the gauge condition implies that the amplitude and the 
wave vector are orthogonal: e°'^kp — 0. 

Linearized theory describes a classical gravitational field whose quan- 
tum description would be a massless spin 2 field that propagates at the 
speed of light. We expect from this that such a field will have only 2 
independent degrees of freedom (helicities in quantum language, polar- 
izations in classical terms). To show this classically we remember that 
hafs is symmetric, so it has 10 independent components, and that Lorentz 
gauge applies 4 independent conditions to these, reducing the freedom to 
6. However, Lorentz gauge doe snot fully fix the coordinates. In fact if we 
perform another infinitesimal coordinate transformation (x'' — )■ x'^ -\- ^^ , 
with ^^^y = 0(h)) and impose □^'^ = 0, we remain in Lorentz gauge. We 
can use this freedom to demand: 

e°" = =^ e^-'fcj = (transverse wave), (1.8) 

e% — 0(traceless wave). (1-9) 

These conditions can only be applied outside a sphere surrounding the 
source. Together they put metric into the transverse-traceless (TT) 
gauge. We will explicitly construct this gauge in chapter 5. 

1.2 Using TT gauge to understand gravitational waves 

The TT gauge leaves only two independent polarizations out of the 

original ten, and it ensures that ha^ — hafj. In order to understand the 
polarization degrees of freedom, let us take the wave to move in z-direction, 
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Figure 1.1. Illustration of two linear polarizations and the associated wave 
amplitude. 
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SO that kz = uj, k^ = uj, k^ = 0, ky = 0; the TT gauge conditions in 
Equation ((TH]) and dO]) lead to e'^" = e^" = and e^^ = -evy. This 
leaves only 2 independent components of the polarization tensor, say e^^ 
and e^y (which we denote by the symbols (B,<E)). 

A wave for which e^^ ~ (pure © polarization) produces a metric of 
the form: 

ds^ = -dt^ + (1 + h+)dx'^ + (1 - h+)d7/ + dz^, (1.10) 

where /i+ = Ae^^ cxp[—iuj{t — z)]. Such a metric produces opposite ef- 
fects on proper distance on the two transverse axes, contracting one while 
expanding the other. If e^^ = we have pure <Si polarization hx which 
can be obtained from the previous case by a simple 45 rotation, as in 
Figure 11.11 Since the wave equation and TT conditions are linear, a gen- 
eral wave will be a linear combination of these two polarization tensors. A 
circular polarization basis would be: 

eR^ -={e++ieyX gl = — =(e+ - iex), (1.11) 

where e_|_, ex are the two linear polarization tensors and e^ and e^ are 
polarizations that rotate in the righthanded and lefthanded directions re- 
spectively. It is important to understand that, for circular polarization, the 
polarization pattern rotates around the central position, but test particles 
themselves rotate only in small circles relative to the central position. 

Now we compute the effects of a wave in the TT gauge on a parti- 
cle at rest in the flat background metric rjap before the passage of the 
gravitational wave. The geodesic equation 

d^xt" dx" dxl^ _ 

dr"^ " dr dr 
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implies in this case: 

^ = -r'oo = -^(2/i»o,o - /loo,.) = 0, (1.12) 

so that the particle does not move. The TT gauge, to first order in 
hai3, represents a coordinate system that is comoving with freely- 
faUing particles. Because hoa = 0, TT-time is proper time on the 
clock of freely-falling particles at rest. 

Tidal forces show the action of the wave independently of the coordi- 
nates. Let us consider the equation of geodesic deviation, which governs 
the separation of two neighboring freely falling test particles A and B. If 
the particles are initially at rest, then as the wave passes it produces an 
oscillating curvature tensor, and the separation ^ of the two particles is: 

^-^'oM'- (1-13) 

To calculate the component R^ojo of Riemann tensor in Equation (I1.13p . 
we can use the metric in the TT gauge, because the Riemann tensor is 
gauge- invariant at linear order (see exercise 4 at the end of this Lecture) . 
So we can replace R^ojo by R^ojo = 5^^"^*j,oo and write: 

This equation, with an initial condition ^L, = constant, describes the oscil- 
lations of B's location as measured in the proper reference frame of A. The 
validity of Equation (|1.14l) is the same as that of the geodesic deviation 
equation: geodesies have to be close to one another, in a neighborhood 
where the change in curvature is small. In this approximation a gravita- 
tional wave is like an extra force, called a tidal force, perturbing the proper 
distance between two test particles. If there are other forces on the par- 
ticles, so that they are not free, then as long as the gravitational field is 
weak, one can just add the tidal forces to the other forces and work as if 
the particle were in special relativity. 

1.3 Interaction of gravitational waves with detectors 

We have shown above that the TT gauge is a particular coordinate system 
in which the polarization tensor of a plane gravitational wave assumes a 
very simple form. This gauge is comoving for freely-falling particles and so 
it is not the locally Minkowskian coordinate system that would be used by 
an experimenter to analyze an experiment. In general relativity one must 
always to be aware of how one's coordinate system is defined. 
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We shall analyze two typical situations: 

• the detector is small compared to the wavelength of the gravitational 
waves it is measuring; and 

• the detector is comparable to or larger than that wavelength. 

In the first case we can use the geodesic deviation equation above to 
represent the wave as a simple extra force on the equipment. Bars detectors 
can always be analyzed in this way. Laser interferometers on the Earth can 
be treated this way too. In these cases a gravitational wave simply produces 
a force to be measured. There is no more to say from the relativity point of 
view. The rest of the detection story is the physics of the detectors. Sadly, 
this is not as simple as the gravitational wave physics! 

In the second case, the geodesic deviation equation is not useful be- 
cause we have to abandon the "local mathematics" of geodesic deviation 
and return to the "global mathematics" of the TT gauge and metric com- 
ponents h'^'^ap- Space-based interferometers like LISA, accurate ranging 
to solar-system spacecraft, and pulsar timing are all in this class. Together 
with ground interferometers, these are beam detectors: they use light 
(or radio waves) to register the waves. 

To study these detectors, it is easiest to remain in the TT gauge and 
to calculate the effect of the waves on the (coordinate) speed of light. Let 
consider, for example, the ® metric from Equation (|1.10p and examine a 
null geodesic moving in the x-direction. The speed along this curve is: 

This is only a coordinate speed, not a contradiction to special relativity. 

To analyze the way in which detectors work, suppose one arm of an 
interferometer lies along the x-direction and the wave, for simplicity, is 
moving in the z-direction with a ® polarization of any waveform h^{t) 
along this axis. (It is a plane wave, so its waveform does not depend on x.) 
Then a photon emitted at time t from the origin reaches the other end, at 
a fixed coordinate position a; = L, at the coordinate time 



tfar^t+ ^l + h+{t{x))dx, (1.16) 

Jo 

where the argument t{x) denotes the fact that one must know the time to 
reach position x in order to calculate the wave field. This implicit equation 
can be solved in linearized theory by using the fact that /i-|_ is small, so we 
can use the first-order solution of Equation (|1.15p above to calculate ft.+ (t) 
to sufficient accuracy. 

To do this we expand the square-root in powers of /i+ , and consider as 
a zero-order solution a photon travelling at speed of light in the a;-direction 
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of a flat space-time. We can set t{x) =t + x. The result is: 

tout^t + L+- h+{t + x)dx. (1.17) 

^ Jo 

In an interferometer, the light is reflected back, so the return trip takes 
1 



'^return t ^ Lj , 



L pL 

h^{t + x)dx + / h^{t + X + L)dx 



(1.18) 



What one monitors is changes in the time taken by a return trip as a 
function of time at the origin. If there were no gravitational wave tretum 
would be constant because L is fixed, so changes indicate a gravitational 
wave. 

The rate of variation of the return time as a function of the start time 
t is 

^^ = l + l[/.+ (t + 2L)-/.+ (t)]. (1.19) 

This depends only on the wave amplitude when the beam leaves and when 
it returns. 

Let us consider now a more realistic geometry than the previous one, 
and in particular suppose that the wave travels at an angle 9 to the z-axis 
in the x-z plane. If we re-do this calculation, allowing the phase of the wave 
to depend on x in an appropriate way, and taking into account the fact that 
^TTxa; jg rcduccd if the wave is not moving in a direction perpendicular to 
X, we find (see exercise 1 at the end of this chapter for the details of the 
calculation) 

^^^ = ^ {(1 - sin^) h%^ (t + 2L) - (1 + sine) hl^it) 

-^2sin0/i^"'[t + L(l-sin6l)]}. (1.20) 

This three-term relation is the starting point for analyzing the response of 
all beam detectors. This is directly what happens in radar ranging or in 
transponding to spacecraft, where a beam in only one direction is used. 
In long-baseline interferometry, one must analyze the second beam as well. 
We shall discuss these cases in turn. 

1.4 Analysis of beam detectors 

1.4.1 Ranging to spacecraft 

Both NASA and ESA perform experiments in which they monitor the re- 
turn time of communication signals with intcr-planetary spacecraft for the 
characteristic effect of gravitational waves. For missions to Jupiter and 
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Saturn, the return times are of order 2-4 x 10'^ s. Any gravitational wave 
event shorter than this wih leave an imprint on the delay time 3 times: once 
when the wave passes the Earth-based transmitter, once when it passes the 
spacecraft, and once when it passes the Earth-based receiver. Searches use 
a form of pattern matching to look for this characteristic imprint. There 
are two dominant sources of noise: propagation-time irregularities caused 
by fluctuations in the solar wind plasma, and timing noise in the clocks 
used to measure the signals. The plasma delays depend on the radio- wave 
frequency, so by using two transmission frequencies one can model and sub- 
tract the plasma noise. Then if one uses the most stable atomic clocks, it 
is possible to achieve sensitivities for h of order 10^^'^. In the future, using 
higher radio frequencies, such experiments may reach 10~^^. No positive 
detections have yet been made, but the chances are not zero. For example, 
if a small black hole fell into the massive black hole in the center of the 
Galaxy, it would produce a signal with a frequency of about 10 mHz and an 
amplitude significantly bigger than 10^^'"'. Rare as this might be, it would 
be a dramatic event to observe. 



1.4.2 Pulsar timing 

Many pulsars, particular old millisecond pulsars, are extraordinarily reg- 
ular clocks, whose random timing irregularities are too small for even the 
best atomic clocks to measure. Other pulsars have weak but observable 
irregularities. Measurements of or even upper limits on any of these tim- 
ing irregularities for single pulsars can be used to set upper limits on any 
background gravitational wave field with periods comparable to or shorter 
than the observing time. Here the 3-term formula is replaced by a simpler 
two-term expression (see exercise 2 at the end of this chapter), because we 
only have a one-way transmission from the pulsar to Earth. Moreover, the 
transit time of a signal to the Earth from the pulsar may be thousands of 
years, so we cannot look for correlations between the two terms in a given 
signal. Instead, the delay time is a combination of the effects of uncorre- 
lated waves at the pulsar when the signal was emitted and at the Earth 
when it is received. 

If one simultaneously observes two or more pulsars, the Earth-based 
part of the delay is correlated between them, and this offers a means of 
actually detecting long-period gravitational waves. Observations require 
timescale of several years in order to achieve the long-period stability of 
pulse arrival times, so this method is suited to looking for strong gravita- 
tional waves with periods of several years. 
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1.4.3 Interferometry 

An interferometer essentially measures changes in the difference in the re- 
turn times along two different arms. It does this by looking for changes 
in the interference pattern formed when the returning light beams are su- 
perimposed on one another. The response of each arm will follow the 
three-term formula in Equation (jl.20p . but with a different value of 9 for 
each arm, depending in a complicate way on the orientation of the arms rel- 
ative to the direction of travel and the polarization of wave. Ground-based 
interferometers are small enough to use the small- L formulas we derived 
earlier. But LISA, the space-based interferometer that will be described 
by Bender at this meeting, is larger than a wavelength of gravitational 
waves for frequencies above 10 mHz, so a detailed analysis of its sensitivity 
requires the full three-term formula. 

Exercises for Lecture 1 

Suggested solutions for these exercises are at the end of the lectures. 
1 (a) Derive the full three-term return equation, reproduced here: 



dt 2 



- {(1 - sine) /i^^ (t + 2L) - (1 + sin 6*) /(,^^(t) 
+2sine'/i^"=[t + L(l-sin6')]} (1.21) 

(b) Show that, in the limit where L is small compared to the wave- 
length of the gravitational wave, the derivative of the return time is the 
derivative of the excess proper distance SL — Lh^ (t) cos^ 9 for small 
L. Make sure you know how to interpret the factor of cos^ 6. 

(c) Examine the limit of the three-term formula when the gravitational 
wave is travelling along the x-axis too (9 — ±-|_j; what happens to 
light going parallel to a gravitational waves 

2 Derive the two-term formula governing the delays induced by gravita- 
tional waves on a signal transmitted only one-way, for example from 
a pulsar to Earth. 

3 A frequently asked question is: if gravitational waves alter the speed 
of light, as we seem to have used here, and if they move the ends of 
interferometer closer and further apart, might these effects not cancel, 
so that there would be no measurable effects on light? Answer this 
question. You may want to examine the calculation above: did we 
make use of the changing distance between the ends, and why or why 
not? 

4 Show that Riemann tensor is gauge-invariant in linearized theory. 



Chapter 2 

Gravitational wave detectors 



Gravitational radiation is a central prediction of general relativity and its 
detection is a key test of the integrity of the theoretical structure of Ein- 
stein's work. But in the long run, its importance as a tool for observational 
astronomy is likely to be even more important. We have excellent obser- 
vational evidence from the Hulse- Taylor binary pulsar system (described 
in Lecture 3) that the predictions of general relativity concerning gravita- 
tional radiation are quantitatively correct. But we have very incomplete 
information from astronomy today about the likely sources of detectable 
radiation. 

The gravitational wave spectrum is completely unexplored, and when- 
ever a new electromagnetic waveband has been opened to astronomy, as- 
tronomers have discovered completely unexpected phenomena. This seems 
just as likely to me to happen again with gravitational waves, especially 
because gravitational waves carry some kinds of information that electro- 
magnetic radiation cannot convey. Gravitational waves are generated by 
bulk motions of masses, and they encode the mass distributions and speeds. 
They are coherent and their low frequencies reflect the dynamical timescales 
of their sources. 

By contrast, electromagnetic waves come from individual electrons ex- 
ecuting complex and partly random motions inside their sources. They 
are incoherent, and individual photons must be interpreted as samples of 
the large statistical ensemble of photons being emitted. Their frequencies 
are determined by microphysics on length scales much smaller than the 
structure of the astronomical system emitting them. From electromagnetic 
observations we can make inferences about this structure only through 
careful modelling of the source. Gravitational waves, by contrast, carry 
information whose connection to the source structure and motion is fairly 
direct. 

A good example is that of massive black holes in galactic nuclei. From 
observations that span the electromagnetic spectrum from radio waves 

10 
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to X-rays, astrophysicists have inferred that black holes of masses up to 
IO^Mq are responsible for quasar emissions and control the jets that power 
the giant radio emission regions. The evidence for the black hole is very 
strong but indirect: no other known object can contain so much mass in 
such a small volume. Gravitational wave observations will tell us about the 
dynamics of the holes themselves, providing unique signatures from which 
they can be identified, measuring their masses and spins directly from their 
vibrational frequencies. The interplay of electromagnetic and gravitational 
observations will enrich many branches of astronomy. 

The history of gravitational wave detection started in the 1960's with 
J. Weber at the University of Maryland. He built the first bar detector: 
it was a massive cylinder of aluminium (^^2x10'^ kg) operating at room 
temperature (300 K) with a resonant frequency of about 1600 Hz. This 
early prototype had a modest sensitivity, around 10^^"^ or 10^^*. 

Despite this poor sensitivity, in the late 1960's Weber announced the 
detection of a population of coincident events between two similar bars 
at a rate far higher than expected from instrumental noise. This news 
stimulated a number of other groups (at Glasgow, Munich, Paris, Rome, 
Bell Labs, Stanford, Rochester, LSU, MIT, Beijing, Tokyo) to build and 
develop bar detectors to check Weber's results. Unfortunately for Weber 
and for the idea that gravitational waves were easy to detect, none of these 
other detectors found anything, even at times when Weber continued to 
find coincidences. Weber's observations remain unexplained even today. 
However the failure to confirm Weber was in a real sense a confirmation of 
general relativity, because theoretical calculations had never predicted that 
reasonable signals would be strong enough to be seen by Weber's bars. 

Weber's announcements have had a mixed effect on gravitational wave 
research. On the one hand, they have created a cloud under which the 
field has labored hard to re-establish its respectability in the eyes of many 
physicists. The legacy of this even today is an extreme cautiousness among 
the major projects, a conservatism that will ensure that the next claim of a 
detection will be ironclad. On the other hand, the stimulus that Weber gave 
to other groups to build detectors has directly led to the present advanced 
state of detector development. 

From 1980 to 1994 groups developed detectors in two different direc- 
tions: 

• Cryogenic bar detectors, developed primarily at Rome/Frascati, 
Stanford, LSU, and Perth (Australia). The best of these detectors 
reach below 10^^^. They are the only detectors operating continuously 
today and and they have performed a number of joint coincidence 
searches, leading to upper limits but no detections. 

• Interferometers, developed at MIT, Garching (where the Munich 
group moved), Glasgow, Caltech, and Tokyo. The typical sensitivity 
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of these prototypes was 10~^^. The first long coincidence observation 
with interferometers was the Glasgow/Garching 100-hour experiment 
in 1989. [2] 

In fact, interferometers had apparently been considered by Weber, but 
at that time the technology was not good enough for this kind of detec- 
tors. Only 10-15 years later, technology had progressed. Lasers, mirror 
coating and polishing techniques, and materials science had advanced far 
enough to allow the first practical interferometers, and it was clear that 
further progress would continue unabated. Soon thereafter several major 
collaborations were formed to build large-scale interferometric detectors: 

• LIGO: Caltech & MIT (NSF) 

• VIRGO: France (CNRS) & Italy (INFN) 

• GEO600: Germany (Max Planck) & UK (PPARC) 

Later other collaborations were formed in Australia (AIGO) and Japan 
(TAMA and JGWO). 

At present there is still considerable effort in building successors to 
Weber's original resonant-mass detector: ultra-cryogenic bars are in oper- 
ation in Frascati and Padova, and they are expected to reach below 10~^°. 
Further, there are proposals for a new generation of spherical or icosahe- 
dral solid- mass detectors from the USA (LSU), Brasil, the Netherlands, 
and Italy. Arrays of smaller bars have been proposed for observing the 
highest frequencies, where neutron star normal modes lie. 

However, the real goal for the near future is to break through the 10^^^ 
level, which is where theory predicts that it is not unreasonable to expect 
gravitational waves of the order of once per year. (See the discussion in 
Lecture 3 below.) The first detectors to reach this level will be the large- 
scale interferometers that are now under construction. They have very long 
arms: LIGO, Hanford (WA) and Livingstone (LA), 4 km; VIRGO: Pisa, 3 
km; GEO600: Hannover, 600 m; TAMA300: Tokyo, 300 m. 

The most spectacular detector in the near future is the space-based 
detector LISA, which has been adopted by ESA (European Space Agency) 
as a Cornerstone mission for the 21** century. The project is now gaining 
a considerable amount of momentum in the USA, and a collaboration be- 
tween ESA and NASA seems likely. This mission could be launched around 
2010. 



2.1 Gravitational wave observables 

We have described earlier how different gravitational wave observables are 
from electromagneric observables. Here are the things that we want to 
measure when we detect gravitational waves: 
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• ^+(i), hx{t), phase{t): the amplitude and polarization of the wave, 
and the phase of polarization, as functions of time. These contain 
most of the information about gravitational wave. 

• 9, (p: the direction on the sky of the source (except for observations of 
a stochastic background). 

From this it is clear that gravitational wave detection is not the same 
as electromagnetic radiation detection. In electromagnetic astronomy one 
almost always rectifies the electromagnetic wave, while we can follow the 
oscillations of the gravitational wave. Essentially in electromagnetism one 
detects the power in the radiation, while for gravitational radiation, as we 
have said before, one detects the wave coherently. 

Let us consider now what we can infer from a detection. If the gravi- 
tational wave has a short duration, of the order of the sampling time of the 
signal stream, then each detector will usually give just a single number, 
which is the amplitude of the wave projected on the detector (a projec- 
tion of the two polarizations /i+ and h^). If the wave lasts more than one 
sampling time, then this information is a function of time. 

If the signal lasts for a sufficiently long time, then both the amplitude 
and the phase of the wave can be affected by the motion of the detector, 
which moves and turns with the motion of Earth. This produces an ampli- 
tude and phase modulation which is not intrinsic to signal. If the signal's 
intrinsic form is understood, then this modulation can be used to determine 
the location of the source. 

We distinguish three distinct kinds of signals, from the point of view 
of observations. 

Bursts have a duration so short that modulation due to detector mo- 
tion is not observable. During the detection, the detector is effectively 
stationary. In this case we need at least 3, and preferably 4, interferome- 
ters to triangulate the positions of bursts on the sky and to find the two 
polarizations h+ and h^. (See discussions in Schutz 1989.) A network of 
detectors is essential to extract all the information in this case. 

Continuous waves by definition last long enough for the motion 
of the detector to induce amplitude and phase modulation. In this case, 
assuming a simple model for the intrinsic signal, we can use the information 
imprinted on the signal (the amplitude modulation and phase modulation) 
to infer the position and polarization amplitude of the source on the sky. A 
single detector, effectively, performs aperture synthesis, finding tne position 
of the source and the amplitude of the wave entirely by itself. However, in 
order to be sure that the signal is not an artefact, it will be important that 
the signal is seen by a second or third detector. 

Stochastic backgrounds can be detected just like noise in a single 
detector. If the detector noise is well understood, this excess noise may be 
detected as a stochastic background. This is closely analogous to the way 
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the original microwave background detection was discovered. 

A more reliable method for detecting stochastic radiation is the cross- 
correlation between two detectors, which experience the same cosmological 
noise but have a different intrinsic noise. Coherent cross-correlation be- 
tween two detectors eliminates much detector noise and works best when 
detectors are closer than a wavelength. 

In general, detection of gravitational waves requires joint observing by 
a network of detectors, both to increase the confidence of the detection and 
to provide accurate information on other physical observables (direction, 
amplitude and so on). Networks can be assembled from interferometers, 
bars, or both. 



2.2 The physics of interferometers 

Interferometric gravitational wave detectors are the most sensitive instru- 
ments, and among the most complex, that have ever been constructed. 
They are remarkable for the range of physics that is important for their 
construction. Interferometer groups work at the forefront of development 
in lasers, mirror polishing and coating, quantum measurement, materials 
science, mechanical isolation, optical system design, and thermal science. 
In this section we shall only be able to take a fairly superficial look at 
one of the most fascinating instrumentation stories of our age. A good 
introduction to interferometer design is Saulson (1994). 

Interferometers use laser light to compare the lengths of two perpen- 
dicular arms. The simplest design, originated by Michelson for his famous 
experiment on the velocity of light, uses light that passes up and down each 
arm once, as in the first panel in Figure 12.11 Imagine such an instrument 
with identical arms defined by mirrors that hang from supports, so they are 
free to move horizontally in response to a gravitational wave. If there is no 
wave, the arms have the same length, and the light from one arms returns 
exactly in phase with that from the other. When the wave arrives, the 
two arms typically respond differently. The arms are no longer the same 
length, and so the light that arrives back at the center from one arm will 
no longer be in phase with that arriving back from the other arm. This will 
produce a shift in the interference fringes between the two beams. This is 
the principle of detection. 

Real detectors are designed to store the light in each arm for longer 
than just one reflection. (See part (b) of Figure [13]) It is optimum to 
store the light for one-half of the period of the gravitational wave, so that 
on each reflection the light gains an added phase shift. Michelson- type 
delay-line interferometers store the light by arranging multiple reflections. 
Fabry-Perot interferometers store the light in cavities in each arm, allowing 
only a small fraction to escape for the interference measurement (part (e) 
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of Figure EJ|) . 

An advantage of interferometers as detectors is that the gravitational- 
wave-induced phase shift of the hght can be made larger simply by making 
the arm-length larger, since gravitational waves act by tidal forces. A 
detector with an arm length I — 4 km responds to a gravitational wave 
with an amplitude of 10~^^ with 

Slgw ^ -hlr^2x 10-^® m (2.1) 

where Slgw is the change in the length of one arm. If the orientation of 
the interferometer is optimum, then the other arm will change by the same 
amount in the opposite direction, so that the interference fringe will shift 
by twice this length. 

If the light path is folded or resonated, as in Panels (b) or (d) of 
Figure 12. 1[ then the effective number of bounces can be traded off against 
overall length to achieve a given desired total path length, or storage time. 
Shorter interferometers with many bounces have a disadvantage, however: 
even though they can achieve the same response as a longer interferometer, 
the extra bounces introduce noise from the mirrors, as discussed below. 
There is therefore a big advantage to long-arm interferometers. 

There are three main sources of noise in interferometers: thermal, shot, 
and vibrational. To understand the way they are controlled, it is important 
to think in frequency space. Observations with ground-based detectors will 
be made in a range from perhaps 10 Hz up to 10 kHz, and initial detectors 
will have a much smaller observing bandwidth within this. Disturbances 
by noise that occur at frequencies outside the observation band can simply 
be filtered out. The goal of noise control is to reduce disturbances in the 
observation band. 

• Thermal noise. Interferometers work at room temperature, and vi- 
brations of the mirrors and of the suspending pendulum can mask 
gravitational waves. To control this noise, scientists take advantage 
of the fact that thermal noise has its maximum amplitude at the fre- 
quency of the vibrational mode, and if the resonance of the mode 
is narrow (a high quality factor Q) then the amplitude at other fre- 
quencies is small. Therefore pendulum suspensions are designed with 
the pendulum frequency at about 1 Hz, well below the observing win- 
dow, and mirror masses are designed to have principal vibration modes 
above 1 kHz, well above the optimum observing frequency for initial 
interferometers. These systems are constructed with high values of Q 
(10^ or more) to reduce the noise in the observing band. Even so, 
thermal noise is typically a dominant noise below 100 or 200 Hz. 

• Shot noise. This is the principal limitation to sensitivity at higher 
frequencies, above 200 to 300 Hz. It arises from the quantization 
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Figure 2.1. Five steps to a gravitational wave interferometer, (a) The simple 
Michelson. Notice that there are two return beams: one goes toward the photode- 
tector and the other toward the laser, (b) Delay line: a Michelson with multiple 
bounces in each arm to enhance the signal, (c) Power recycling. The extra mirror 
recycles the light that goes towards the laser, which would otherwise be wasted, 
(d) Signal recycling. The mirror in front of the photodetector recycles only the 
signal sidebands, provided that in the absence of a signal no light goes to the 
photodetector. (e) Fabry-Perot interferometer. The delay lines are converted to 
cavities with partially silvered interior mirrors. 
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of photons. When photons form interference fringes, they arrive at 
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random times and make random fluctuations in the light intensity that 
can look like a gravitational wave signal; the more photons one uses, 
the smoother will be the interference fringe. We can easily calculate 
this intrinsic noise. If N is the number of photons emitted by the laser 
during our measurement, then as a random process the fluctuation 
number 5N is proportional to the square root of N. If we are using 
light with a wavelength A (for example infrared light with A ~ 1 /zm) 
one can expect to measure lengths to an accuracy of 

Olshot 



2ttVN 

To measure a gravitational wave at a frequency /, one has to make at 
least 2/ measurements per second, so one can accumulate photons for 
a time i. If P is the light power, one has 



N 



he J_ 
A • 2/ 



It is easy to work out from this that, for Slshot to be equal to Slgw in 
Equation (j2.ip , one needs light power of about 600 kW. No continuous 
laser could provide this much light to an interferometer. 
The key to reaching such power levels inside the arms of a detector is 
a technique called power recycling (see Saulson 1994), first proposed 
by Drever and independently by Schilling. Normally, interferometers 
work on a "dark fringe", that is they are arranged so that the light 
reaching the photodetector is zero if there is no gravitational wave. 
Then, as shown in part (a) of Figure [231 the whole of the input light 
must emerge from the interferometer travelling towards the laser. If 
one places another mirror, correct positioned, between the laser and 
the beam splitter (part (c) of the figure), it will reflect this wasted light 
back into the interferometer in such a way that it adds coherently in 
phase with light emerging from laser. In this way, light can be recycled 
and the required power levels in the arms achieved. 
Of course, there will be a maximum recycling gain, which is set by mir- 
ror losses. Light power builds up until the laser merely resupplies the 
losses at the mirrors, due to scattering and absorption. The maximum 
power gain is 

P ' 



l-i?2 



where 1 — i?^ is the total loss summed over all the optical surfaces. For 
the very high-quality mirrors used in these projects, 1 — R^ ~ 10~^. 
This reduces the power requirement for the laser by the same factor, 
down to about 6 W. This is attainable with modern laser technology. 
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• Ground vibration and mechanical vibrations are another source 
of noise that must be screened out. Typical seismic vibration spec- 
tra fall sharply with frequency, so this is a problem primarily below 
100 Hz. Pendulum suspensions are excellent mechanical filters above 
the pendulum frequency: it is a familiar elementary-physics demon- 
stration that one can wiggle the suspension point of a pendulum vig- 
orously at a high frequency and the pendulum itself remains undis- 
turbed. Suspension designs typically involve multiple pendula, each 
with a frequency around 1 Hz. These provide very fast roll-off of 
the noise above 1 Hz. Interferometer spectra normally show a steep 
low-frequency noise "wall" : this is the expected vibrational noise am- 
plitude. 

In addition, there are noise sources that are not dominant in the 
present interferometers but will become important as sensitivity increases. 

• Quantum effects: uncertainty principle noise. Shot noise is 
a quantum noise, but in addition there are other effects similar to 
those that bar detectors face, as described below: zero-point vibrations 
of suspensions and mirror surfaces, and back-action of light pressure 
fluctuations on the mirrors. These are small compared to present 
operating limits of detectors, but they may become important in 5 
years or so. Practical schemes to reduced this noise have already been 
demonstrated in principle, but they need to be improved considerably. 
This is the subject of considerable theoretical work at the moment. 

• Gravity gradient noise. Gravitational wave detectors respond to 
any changes in the gradients (tidal forces) of the local gravitational 
field, not just those carried by waves. The environment always con- 
tains changes in the Newtonian fields of nearby objects. Besides ob- 
vious ones, like people, there are changes caused by density waves 
in ground vibrations, atmospheric pressure changes, and many other 
disturbances. Below about 1 Hz, these gravity gradient changes will 
be stronger than waves expected from astronomical objects, and they 
make it impossible to do observing at low frequencies from Earth. This 
is the reason that scientists have proposed the LISA mission, discussed 
below. Above 1 Hz, this noise does not affect the sensitivity of present 
detectors, but in 10 years this could become a limiting factor. 

Besides these noise sources, which are predictable and therefore can be 
controlled by detector design, it is possible that there will be unexpected or 
unpredicted noise sources. Interferometers will be instrumented with many 
kinds of environmental monitors, but there may occasionally be noise that 
is impossible to identify. For this reason, short bursts of gravitational 
radiation must be identified at two or more separated facilities. Even if 
detector noise is not at all understood, it is relatively easy to estimate from 
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the observed noise profile of the individual detectors what the chances are 
of a coincident noise event between two detectors. 



2.2.1 New Interferometers and their capabilities 

Interferometers work over a broad bandwidth and they do not have any 
natural resonance in their observing band. They are ideal for detecting 
bursts, since one can perform pattern-matching over the whole bandwidth 
and detect such signals optimally. They are also ideal for searching for 
unknown continuous signals, such as surveying the sky for neutron stars. 
And in observations of stochastic signals by cross-correlating two detectors, 
they can give information about the spectrum of the signal. 

If an interferometer wants to study a signal with a known frequency, 
such as a known pulsars, then there is another optical technique available to 
enhance its sensitivity in a narrow bandwidth, at the expense of sensitivity 
outside that band. This is called signal recycling[Z . In this technique, 
a further mirror is placed in front of the photo-detector, where the signal 
emerges from the interferometer (see part (d) of Figure l2.1l) . If the mirror is 
chosen correctly, it will build up the signal, but only in a certain bandwidth. 
This modifies the shot noise in the detector, but not other noise sources. So 
it can improve sensitivity only at the higher frequencies where shot noise 
is the limiting factor. 

Four major interferometer projects are now under construction, and 
they could begin acquiring good data in the period between 2000 and 2003. 
They will all operate initially with a sensitivity approaching 10~^^ over a 
bandwidth between 50 and 1000 Hz. Early detections are by no means cer- 
tain, but recent work has made prospects look better for an early detection 
than when these detectors were funded. 

TAMA300 4 (Japan) is located in Tokyo, and its arm length is 300 m. 
It began taking data without power recycling in 1999, but its sensitivity is 
not yet near 10~^^. Following improvements, especially power recycling, it 
should get to within a factor of 10 of this goal. But it is not planned as an 
observing instrument: it is a prototype for a kilometer-scale interferometer 
in Japan, currently called JGWO. By 2005 this may be operating, possibly 
with cryogenically cooled mirrors. 

GEOBOOfS" (Germany & Britain) is located near Hannover (Ger- 
many). Its arm length is 600 m and the target date for first good data 
is now the end of 2001. Unlike TAMA, GEO600 is designed as a leading- 
edge-technology detector, where high-performance suspensions and optical 
tricks like signal recycling can be developed and applied. Although it has 
a short baseline, it will have a similar sensitivity to the larger LIGO and 
VIRGO detectors at first. At a later stage, LIGO and VIRGO will incor- 
porate the advanced methods developed in GEO, and at that point they 
will advance in sensitivity, leaving GEO behind. 
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Figure 2.2. TAMA300 sensitivity as a function of frequency. The vertical axis 
is the 1(T noise level, measured in strain per root Hz. To get a limit on the 
gravitational wave amplitude h, one must multiply the height of the curve by the 
square root of the bandwidth of the signal. This takes account of the fact that the 
noise power at different frequencies is independent, so the power is proportional 
to bandwidth. The noise amplitude is therefore proportional to the square root 
of the bandwidth. 
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Figure 2.3. GEO600 noise curves. As for the TAMA curve, these are calibrated 
in strain per root Hz. The figure on the left shows GEO's wideband configuration; 
that on the right shows a possible narrowband operating mode. 
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As we can see from Figure [^751 that the sensitivity of GEO600 depends 
on its bandwidth, which in its turn depends on the signal recycling factor. 
GEO600 can change its observing bandwidth in response to observing goals. 
By choosing low or high reflectivity for the signal recycling mirror, scientists 
can make GEO600 wide-band or narrow-band, respectively. The center 
frequency of the observing band (in the right panel of Figure 12.31 it is ^ 
600 Hz) can be tuned to any desired frequency by shifting the position of 
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Figure 2.4. Noise curves of the initial LIGO (left) and VIRGO (right) detectors. 
The VIRGO curve is in strain per root Hz, as the GEO curves earlier. The LIGO 
curve is calibrated in meters per root Hz, so to convert to a limit on h one 
multiplies by the square root of the bandwidth and divides by the length of the 
detector arm, 4000 m. 




\l ^ 11,1 



the signal recycling mirror, thus changing the resonance frequency of the 
signal recycling cavity. This feature could be useful when interferometers 
work with bars or when performing wide-band surveys. 

LIGO [6 (USA) is building two detectors of arm length 4 km. One is 
located in Hanford WA and the other in Livingstone LA. The target date for 
observing is mid-2002. The two detectors are placed so that their antenna 
patterns overlap as much as possible and yet they are far enough apart 
that there will be a measurable time-delay in most coincident bursts of 
gravitational radiation. This delay will give some directional information. 
The Hanford detector also contains a half-length interferometer to assist in 
coincidence searches. The two LIGO detectors are the best placed for doing 
cross-correlation for a random background of gravitational waves. LIGO's 
expected initial noise curve is shown in Figure 12.41 These detectors have 
been constructed to have a long lifetime. With such long arms they can 
benefit from upgrades in laser power and mirror quality. LIGO has defined 
an upgrade goal called LIGO II, which it hopes to reach by 2007, which 
will observe at 10^^^ or better over a bandwidth from 10 Hz up to 1 kHz. 

VIRGO [7 (Italy & France) is building a 3 km detector near Pisa. 
Its target date for good data is 2003. Its expected initial noise curve is 
shown in Figure 12.41 Like LIGO, it can eventually be pushed to much 
higher sensitivities with more powerful lasers and other optical enhance- 
ments. VIRGO specializes in sophisticated suspensions, and the control of 
vibrational noise. Its goal is to observe at the lowest possible frequencies 
from the ground, at least partly to be able to examine as many pulsars and 
other neutron stars as possible. 
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2.3 The physics of Resonant Mass Detectors 

The principle of operation of bar detectors is to use the gravitational tidal 
force of the wave to stretch a massive cylinder along its axis, and then to 
measure the elastic vibrations of the cylinder. 

Let us suppose we have a typical bar with length L ^ 1 ra. (In the 
future, spheres may go up to 3 m.) Depending on the length of the bar 
and its material, the resonant frequency will be / ^^ 500 Hz to 1.5 kHz and 
mass M ^ 1000 kg. A short burst gravitational wave h will make the bar 
vibrate with an amplitude 



Slaw ^hl^ 10^21 



•'gw 



m. 



Unlike the interferometers, whose response is simply given by this equation, 
the bars respond in a complicated way depending on all their internal 
forces. But if the duration of the wave is short, the amplitude will be of 
the same order as that given here. If the wave has long duration and is 
near the bar's resonant frequency, then the signal can build up to much 
larger amplitudes. Normally, bar detector searches have been targeted at 
short-duration signals. 

The main sources of noise that compete with this very small amplitude 



are: 



Thermal noise. This is the most serious source of noise. Interfer- 
ometers can live with room-temperature thermal noise because their 
larger size makes their response to a gravitational wave larger, and 
because they observe at frequencies far from the resonant frequency, 
where the noise amplitude is largest. But bars observe at the resonant 
frequency and have a very short length, so they must reduce thermal 
noise by going to low temperatures. The best ultra-cryogenic bars 
today operate at about T = 100 mK, where the r.m.s. amplitude of 
vibration is found by setting the kinetic energy of the normal mode, 
M((5/)^/2, equal to fcT/2, the equipartition thermal energy of a single 
degree of freedom. This gives then 



(si^)^ = ( — r - 6 • 10"^^ 



m 



This is far larger than the gravitational wave amplitude. In order to 
detect gravitational waves against this noise, bars are constructed to 
have a very high Q, of order 10^ or better. 

The reason that bars need a high Q is different from the reason that 
interferometers also strive for high-Q systems. To see how Q helps 
bars, recall that Q is defined as Q = f ■ t where / is the resonant 
frequency of the mode and r is the decay time of the oscillations. If Q 
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is large, then the decay time is long. If the decay time is long, then the 
amplitude of oscillation changes very slowly in thermal equilibrium. 
Essentially, the bar's mode of vibration changes its amplitude by a 
random walk with very small steps, taking a time H ~ 1000 s to change 
by the full amount. On the other hand, a gravitational wave burst will 
cause an amplitude change in a time of order 1 ms, during which 
the thermal noise will have random walked to an expected amplitude 

change that is Q^ — (-p^) ^ times smaller. In this case 



1 / UJ^ \ 2 
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Thus, thermal noise only affects a measurement to the extent that it 
changes the amplitude of vibration during the time of the gravitational 
wave burst. This change is similar to that produced by a gravitational 
wave of amplitude 6 x 10^^^. It follows that, if thermal noise were 
the only noise source, bars would be operating at around lO^^" today. 
Bar groups expect in fact to reach this level during the next few years, 
as they reduce the other competing sources of noise. Notice that the 
effect of thermal noise has nothing to do with the frequency of the 
disturbance, so it is not the reason that bars observe near their reso- 
nant frequency. In fact, both thermal impulses and gravitational wave 
forces are mechanical forces on the bar, and the ratio of their induced 
vibrations is the same at all frequencies for a given applied impulsive 
force. 

Sensor noise. Because the oscillations of the bar are very small, bars 
require a transducer to convert the mechanical energy of vibration into 
electrical energy, and an amplifier that increases the electrical signal to 
a level where it can be recorded. If the amplifier were perfect, then the 
detector would in fact be broad-band: it would amplify the smaller off- 
resonant responses just as well as the on-resonance ones. Conversely, 
real bars are narrow-bad because of sensor noise, not because of their 
mechanical resonance. 

Unfortunately sensing is not perfect: amplifiers introduce noise and 
this makes small amplitudes harder to measure. The amplitudes of 
vibration are largest in the resonance band near the resonant frequency 
/o, so amplifier noise limits the detector sensitivity to frequencies near 
/q. Now, the signal (a typical gravitational wave burst) has a duration 
time T^ ^ 1 ms, so the amplifier's bandwidth should be at least 1/tiu 
in order for it to be able to record a signal every r^. In other words, 
bars require amplifiers with very small noise in a large bandwidth (^ 
1000 Hz) near /o (note that this band is much larger than f/Q). 
Today typical band widths of realizable amplifiers are 1 Hz, but in the 
very near future it is hoped to extend these to 10 Hz, and eventually 
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to 100 Hz. 

Quantum limit. According to the Heisenberg uncertainty principle, 

the zero-point vibrations of a bar with a frequency of 1 kHz have r.m.s. 

amplitude 

/^/2\^ f ^^ y ^ 4 . 10-21 m 

\ /quant \2ttM f J 

This is bigger than the expected signal, and comparable to the thermal 
limit over 1 ms. It represents the accuracy with which one can measure 
the amplitude of vibration of the bar. So as soon as current detectors 
improve their thermal limits, they will run into the quantum limit, 
which must be overcome before a signal at IQ-^^ can be seen with 
such a detector. One way to overcome this limit is by increasing the 
size of the detector and even by making it spherical. This increases its 
mass dramatically, pushing the quantum limit down below lO-^^. 
Another way around the quantum limit is to avoid measuring 51 ^ but 
instead to measure other observables. After all, the goal is to infer 
the gravitational wave amplitude, not to measure the state of vibra- 
tion of the bar. It is possible to define a pair of conjugate observables 
that have the property that one of them can be measured arbitrarily 
accurately repeatedly, so that the resulting inaccuracy of knowing the 
conjugate variable's value does not disturb the first variable's value. 
Then, if the first variable responds to the gravitational wave, the grav- 
itational wave may be measured accurately, even though the full state 
of the bar is poorly known. This method is called "hack reaction eva- 
sion". The theory was developed in a classic paper by Caves et al[5]. 
However, no viable schemes to do this have been demonstrated for bar 
detectors so far. 



2.3.1 New Bar Detectors and their capabilities 

Resonant-mass detectors are limited by properties of materials and, as 
we have just explained, they have their best sensitivity in a narrow band 
around their resonant frequency. However they can usefully explore higher 
frequencies (above 500 Hz), where the interferometer noise curves are rising 
(see earlier figures). 

From the beginning, bars were designed to detect bursts. If the burst 
radiation carries significant energy in the bar's bandwidth, then the bar 
can do well. Standard assumptions about gravitational collapse suggest 
a signal with a broad spectrum to 1 kHz or more, so that most of the 
sensitive bars today would be suited to observe such a signal. Binary 
coalescence has a spectrum that peaks at low frequencies, so bars are not 
partiularly well-suited for such signals. On the other hand, neutron-star 
and stellar-mass black-hole normal modes range in frequency from about 
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1 kHz up to 10 kHz, so suitably designed bars could in principle go after 
these interesting signals. 

A bar gets all of its sensitivity in a relatively narrow bandwidth, so if 
a bar and an interferometer can both barely detect a burst of amplitude 
10~^°, then the bar has much greater sensitivity than the interferometer 
in its narrow band, and much worse at other frequencies. This has led 
recently to interest in bars as detectors of continuous signals. If the 
signal frequency is in the observing band of the bar, it can do very well 
compared to interferometers. Signals from millisecond pulsars and possible 
signals from X-ray binaries are suitable if they have the right frequency. 
But most known pulsars will radiate at frequencies rather low compared to 
the operating frequencies of present-day bars. 

The excellent sensitivity of bars in their narrow bandwidth also suits 
them to detecting stochastic signals. Cross-correlations of two bars or of 
bars with interferometers can be better than searches with first-generation 
interferometers. [9 One gets no spectral information, of course, and in the 
long run expected improvements in interferometers will overtake bars in 
this regard. 

Today's best bar detectors are orders of magnitude more sensitive 
than the original Weber bar. Two ultra-cryogenic bars have been built 
and are operating at thermodynamic temperatures below 100 mK: NAU- 
TILUS ^lOj at Frascati, near Rome, and AURIGALll in Legnaro. With 
a mass of several tons, these may be the coldest massive objects ever seen 
anywhere in the universe. These are expected soon to reach a sensitivity 
of 10"^" near 1 kHz. Already they are performing coincidence experiments 
with bars at around 4 K at Perth, Australia, and at LSU. 

Proposals exist in the Netherlands, Brazil, Italy, and the USA for 
spherical or icosahedral detectors. (See links from reference [lOJ.) 
These detectors have more mass, so they could reach 10~^^ near 1 kHz. 
Because of their shape, they have omni-directional antenna patterns; if 
they are instrumented so that all 5 independent fundamental quadrupolar 
modes of vibration can be monitored, they can do all-sky observing and 
determine directions as well as verify detections using coincidences between 
modes of the same antenna. 



2.4 A detector in space 

As we have noted earlier, gravitational waves from astronomical objects at 
frequencies below 1 Hz are obscured by Earth-based gravity-gradient noise. 
Detectors must go into space to observe in this very interesting frequency 
range. 

The LISA [12] mission is likely to be the first such mission to fly. LISA 
will be a triangular array of spacecraft, with arm lengths of 5 x 10^ km, 



A detector in space 26 

orbiting the Sun in the Earth's orbit, about 20° behind the Earth. The 
spacecraft wih be in a plane inchned to the echptic by 60°. The three arms 
can be combined in different ways to form two independent interferometers. 
During the mission the configuration of spacecraft rotates in its plane, and 
the plane rotates as well, so that LISA's antenna pattern sweeps the sky. 

LISA has been named a Cornerstone mission of the European Space 
Agency (ESA), and NASA has recently formed its own team to study the 
same mission, with a view toward a collaboration with ESA. LISA will be 
sensitive in a range from 0.3 mHz to about 0.1 Hz, and it will be able to 
detect known binary star systems in the Galaxy and binary coalescneces 
of supermassive black holes anywhere they occur in the Universe. A joint 
ESA-NASA project looks very likely, aiming at a launch around 2010. A 
technology demonstration mission might be launched in 2005 or 2006. 

LISA'S technology is fascinating. We can only allude to the most 
interesting parts of the mission here. A full description can be found in the 
pre-Phase A study document. 1"5] The most innovative aspect of the mission 
is drag-free control. In order to guarantee that the interferometry is not 
disturbed by external forces, such as fluctuations in solar radiation pressure, 
the mirror that is the reference point for the interferometry is on a free mass 
inside the spacecraft. The spacecraft acts as an active shield, sensing the 
position of the free mass, firing jets to counteract external forces on itself 
and ensure that it does not disturb the free mass. The jets themselves are 
remarkable, in that they must be very weak compared to most spacecraft's 
control jets, and they must be capable of very precise control. They will 
work by expelling streams of ions, accelerated and controlled by a high- 
voltage electric field. Fuel for these jets is not a problem: one gram will be 
enough for a mission lifetime of 10 years! 

LISA interferometry is not done with refiection from mirrors. When a 
laser beam reaches one spacecraft from the other, it is too weak to reflect: 
the sending spacecraft would only get the occasional photon! Instead, the 
incoming light is sensed, and an on-board laser is slaved to it, returning an 
amplified beam with the same phase and frequency as the incoming one. No 
space mission has yet implemented this kind of laser-transponding. The 
LISA team had to ensure that there was enough information in all the 
signals to compensate for inevitable frequency fluctuations among all six 
on-board lasers. 

A further serious problem that the LISA team had to solve was how 
to compensate for the relative motions of the spacecraft. The laser signals 
converging on a single spacecraft from the other two corners will be Doppler 
shifted so that their fringes change at MHz frequencies. This has to be 
sensed on board and removed from the signal that is sent back to Earth, 
which can only be sampled a few tens of times per second. 

When LISA flies it will, on a technical as well as a scientific level, be 
a worthy counterpart to its Earth-based interferometer cousins! 
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Figure 2.5. LISA sensitivity to binary systems in the Galaxy (top) and to mas- 
sive black hole coalescences (bottom). The top figure is calibrated in the intrinsic 
amplitude of the signal, and the noise curve shows the detection threshold {5a) 
for a 1-year observation. It also shows the confusion limit due to unresolved 
binary systems. The bottom panel shows the effective amplitude of signals from 
coalescences of massive black holes. Since some such events last less than 1 year, 
what is shown is the expected signal-to-noise ratio of the observation. 
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2.4.1 LISA'S capabilities 



In the low-frequency LISA window, most sources will be relatively long- 
lived, at least a few months. During an observation, LISA will rotate and 
change its velocity by a significant amount. This will induce Doppler shifts 
into the signals, and modulate their amplitudes, so that LISA should be 
able to infer the position, polarization, and amplitude of sources entirely 
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from its own observations. Below about 1 mHz, this information weakens, 
because the wavelength of the radiation becomes comparable to or greater 
than the radius of LISA's orbit. The amplitude modulation is the only 
directional information in this frequency range. 

2.5 Gravitational and electromagnetic waves compared 
and contrasted 

To conclude this Lecture is useful to discuss the most important differences 
and similarities between gravitational waves and electromagnetic ones. We 
do this in the form of a table. 



Electromagnetism 



General Relativity 



Two signs of charges - large bodies One sign of mass - gravity accumu- 



usually neutral - waves usually emit- 
ted by single particles, often incoher- 
ently - waves carry "local" informa- 
tion. 



lates - waves emitted more strongly 
by larger body - waves carry "global" 
information. 



A genuine physical force, acting dif- 
ferently on different bodies. Detected 
by measuring accelerations. 



Equivalence principle: gravity affects 
all bodies in the same way. Rep- 
resented as a space-time curvature 
rather than a force. Detected only 
by tidal forces - differential accelera- 
tions. 



Ma:xweirs equations are linear. Phys- 
ical field is F^i, (E and B). Gauge 
field is vector potential A. 



Einstein's equations are nonlinear. 
Physical field is Riemann curvature 
tensor Rfivap. Gauge fields are met- 
ric Qfj^v and connection F"^. Gauge 
transformations are coordinate trans- 
formations. 



Source is charge-current density J^. 
Charge creates electric field, current 
magnetic field. 



Source is stress-energy tensor T^^. 
Mass creates a Newtonian-like field, 
momentum as gravito-magnetic ef- 
fects. Stress creates field too. 



Moderately strong force on the 
atomic scale: ip^ = io39. 



Weaker than "weak" interaction. 



Wave generation 

d'^dpA^ = AiveoJix in 
gauge (Lorentz gauge). 



for Ap: 

a convenient 



Wave generation for h^^ 
convenient gauge. 



' 3mi' '7m'^- 



SttT, 



flL- 



in a 
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Propagate at the speed of light, am- Propagate at the speed of hght, am- 
phtude falls as 1/r. plitude falls as 1/r. 



Conservation of charge => radiation 
by low-velocity charges is dominated 
by dipole component. 



Conservation of mass and momentum 
=> radiation by low-velocity masses 
is dominated by quadrupole compo- 
nent. 



Simple detector: oscillating charge. 
Action is along a line, transverse to 
the directions of propagation. Spin 
s = 1 and two states of linear polari- 
sation that are inclined each other at 
an angle of 90° . 



Simple detector: distorted ring of 
masses. Action is elliptic in a plane 
transverse to the direction of propa- 
gation. Spin s = 2 and two states 
of linear polarisation that are inclined 
each other at an angle of 45° . Equiva- 
lence principle ^ action depends only 



Strength of force ; 
refract easily. 



■ waves scatter and 



Weakness of gravity => waves prop- 
agate almost undisturbed and trans- 
fer energy very weakly. Dimension- 
less amplitude h is small. 



Local energy and flux well-defined: 
Poynting vector etc. 



Equivalence principle => local energy 
density cannot be defined exactly. 
Only global energy balance is exact. 



Multipole expansion in slow-motion 
limit is straightforward, radiation re- 
action well-defined. 



Multipole expansion different if fields 
are weak or strong. For quasi- 
Newtonian case fields are weak, and 
the resulting post-Newtonian expan- 
sion is delicate. Radiation reaction is 
still not fully understood. 



Exact solutions, containing waves, 
are available and can guide the con- 
struction of approximation methods 
for more complicated situations. 



Fully realistic exact solutions for dy- 
namical situations of physical interest 
are not available. Extensive reliance 
on approximation methods. 



Chapter 3 

Astrophysics of gravitational wave 
sources 



There are a large number of possible gravitational wave sources in the 
observable wave band, which spans 8 orders of magnitude in frequency: 
from lO"** Hz (lower bound of current space-based detector designs) to 
10^ Hz (frequency limit of likely ground-based detectors). Some of these 
sources are highly relativistic and not too massive, especially above 10 Hz: 
a black hole of mass 1000 M© has a characteristic frequency of 10 Hz, 
and larger holes have lower frequencies in inverse proportion to the mass. 
Neutron stars have even higher characteristic frequencies. Other systems 
are well-described by Newtonian dynamics, such as binary orbits. 

For nearly-Newtonian sources the post-Newtonian approximation (see 
Lecture 5) provides a good framework for calculating gravitational waves. 
More relativistic systems, and unusual sources like the early universe, re- 
quire more sophisticated approaches (see Lecture 6). 

3.1 Sources detectable from ground and from space 

3.1.1 Supernovae and gravitational collapse 

The longest-expected and still probably the least-understood source, grav- 
itational collapse is one of the most violent events know to astronomy. Yet, 
because we have little direct information about the deep interior, we cannot 
make reliable predictions about the gravitational radiation from it. 

Supernovae are triggered by the gravitational collapse of the interior 
degenerate core of an evolved star. According to current theory the result 
should be a neutron star or black hole. The collapse releases an enormous 
amount of energy, about O.ISMqc^, most of which is carried away by neu- 
trinos; an uncertain fraction is converted into gravitational waves. One 
mechanism for producing this radiation could be dynamical instabilities in 
the rapidly rotating core before it becomes a neutron star. Another likely 
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source of radiation is the r-mode instability (see Lecture 6). This could 
release ^ O.IMqc^ in radiation every time a neutron star is formed. 

However, both kinds of mechanisms are difficult to model. The prob- 
lem with gravitational collapse is that perfectly spherical motions do not 
emit gravitational waves, and it is still not possible to estimate in a reli- 
able way the amount of asymmetry in gravitational collapse. Even modern 
computers are not able to perform realistic simulations of gravitational 
collapse in 3D, including all the important nuclear reactions and neutrino- 
and photon-transport. Similarly, it is hard to model the r-mode instability 
because its evolution depends on nonlinear hydrodynamics and on poorly 
known physics, such as the cooling and viscosity of neutron stars. 

An alternative approach is to use general energy considerations. If 
for example we assume that 1% of the available energy is converted into 
gravitational radiation, then, from formulas we will derive in the next Lec- 
ture, the amplitude h would be large enough to be detected by the first 
ground-based interferometers (LIGO/GEO600/VIRGO) at the distance of 
Virgo Cluster (18 Mpc) if the emission centers at 300 Hz. Moreover, bar 
and spherical-mass detectors with an effective sensitivity of 10^^^ and the 
right resonant frequency could see these signals as well. 

The uncertainties in our predictions have a positive aspect: it is clear 
that if we can detect radiation from supcrnovae, we will learn much that we 
don't know about the end stages of stellar evolution and about neutron-star 
physics. 



3.1.2 Binary stars 

Binary systems have given us our best proof of the reliability of general 
relativity for gravitational waves. The most famous example of such sys- 
tems is the binary pulsar PSR1916-f 16, discovered by Hulse and Taylor in 
1974; they were awarded the Nobel Prize for this discovery in 1993. From 
the observations of the modulation of the pulse period as the stars move in 
their orbits, one knows many important parameters of this system (orbital 
period, eccentricity, masses of the two stars, etc), and the data also show 
directly the decrease of the orbital period due to the emission of gravita- 
tional radiation. The observed value is 2.4 • lO"'^^ s/s. Post-Newtonian 
theory allows one to predict this from the other measured parameters of 
the system, without any free parameters (see Lecture 6); the prediction is 
2.38 • 10~^^, in agreement within the measurement errors. 

Unfortunately the radiation from the Hulse- Taylor system will be too 
weak and of too low frequency to be detectable by LISA. 
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3.1.3 Chirping binary systems 

If a binary gives off enough, energy for its orbit to shrink by an observable 
amount during an observation, it is said to chirp: as the orbit shrinks, the 
frequency and amplitude go up. LISA will see a few chirping binaries. If a 
binary system is compact enough to radiate above 10^^ Hz, it will always 
chirp within one year, provided its components have a mass above about 
1 Mq. If they are above about 10'^ Mq, the binary will go all the way to 
coalescence within the one-year observation. 

Chirping binary systems are more easily detectable than gravitational 
collapse events because one can model with great accuracy the gravitational 
wave-form during the inspiral phase. There will be radiation, possibly with 
considerable energy, during the poorly understood plunge phase (when the 
objects reach the last stable orbit and fall rapidly towards one another) 
and during the merger event, but the detectability of such systems rests on 
tracking their orbital emissions. 

The major uncertainty about this kind of source is the event rate. 
Current pulsar observations suggest that there will be ~ 1 coalescence per 
year of a Hulse- Taylor binary out to about 200 Mpc. This is a lower limit 
on the event rate, since it comes from systems we actually observe. It 
is possible that there are other kinds of binaries that we have no direct 
knowledge of, which will boost the event rate. 

Theoretical modeling of binary populations gives a wide spectrum of 
mutually inconsistent predictions. Some authors [H] suggest that there may 
be a large population that escapes pulsar surveys but brings the nearest 
neutron star coalescence in one year as far as 30 Mpc, only slightly fur- 
ther than the Virgo cluster; but other models[15j put the rate near to the 
observational limit. 

The most exciting motivation for detecting coalescing binaries is that 
they could be associated with gamma-ray bursts. The event rates are 
consistent, and neutron stars are able to provide the required energy. If 
gamma-bursts are associated with neutron-star coalescence, then observa- 
tions of coalescence radiation should be followed within a second or so by 
a strong gamma-ray burst. 

LISA will see a few chirping binaries in the Galaxy, but the sensitivity 
of the first generation of ground based detectors is likely to be too poor to 
see many such events (see Table 13.11) . 

A certain fraction of such system could contain black holes instead 
of neutron stars. In fact black holes should be over-represented in binary 
systems (relative to their birth rate) because their formation is much less 
likely to disrupt a binary system (there is much less mass lost) than the 
formation of a neutron star would be. Pulsar observations have not yet 
turned up a black-hole/neutron-star system, and of course one does not 
expect to see binary black holes electromagnetically. So we can only make 
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Table 3.1. The range for detecting a 2 x 1.4 M© NS binary coalescence. The 
threshold for detection is taken to be 5a. The binary and detector orientations 
are assumed optimum. The average S/N ratio for randomly oriented systems is 
reduced from the optimum by 1/VE. 

Detector: TAMA300 GEO 600 LIGO I VIRGO LIGO II 

Range 3 Mpc 14 Mpc 30 Mpc 36 Mpc 500 Mpc 

(S/N=5) 



theoretical estimates, and there are big uncertainties. 

Some evolution calculations [13] suggest that the coalescence rate of 
BH-BH systems may be of the same order as the NS-NS rate. Other 
models[TS] suggest it could even be zero, because stellar-wind mass loss 
(significant in very massive stars) could drive the stars far apart before the 
second BH forms, leading to coalescence times longer than the age of the 
Universe. A recent proposal identifies globular clusters as "factories" for 
binary black holes, forming binaries by 3-body collisions and then expelling 
them. [TB] Gamma-ray bursts may also come from black- hole/neutron-star 
coalescences. If the more optimistic event rates are correct, then black-hole 
coalescences may be among the first sources detected by ground-based de- 
tectors (Table [ 



Table 3.2. The range for detecting a 10 M0 black-hole binary coalescence. 
Conventions as in Table ISTTI 



Detector: GEO 600 LIGO I VIRGO LIGO II 

Range (S/N=5) 75 Mpc 160 Mpc 190 Mpc 2.6 Gpc 



3.1.4 Pulsars and other spinning neutron stars 

There are a number of ways in which a spinning neutron star may give off 
a continuous stream of gravitational waves. They will be weak, so they will 
require long continuous observation times, up to many months. Here are 
some possible emission mechanisms for neutron stars. 

The r-modes. Neutron stars are born hot and probably rapidly ro- 
tating. Before they cool (during their first year) they have a family of 
unstable normal modes, the r-modes. These modes are excited to insta- 
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bility by the emission of gravitational radiation, as predicted originally by 
Andersson [17] . They are particularly interesting theoretically because the 
radiation is gravito-magnetic, generated by mass currents rather than mass 
asymmetries. We will study the theory of this radiation in Lecture 5. In 
Lecture 6 we will discuss how the emission of this radiation excites the 
instability (the CFS instability mechanism). 

Being unstable, young neutron stars will presumably radiate away 
enough angular momentum to reduce their spin and become stable. This 
could lower the spin of a neutron star to ~ 100 Hz within one year after its 
formation [18]. The energy emitted in this way should be a good fraction 
of the star's binding energy, so in principle this radiation could be detected 
from the Virgo Cluster by LIGO II, provided matched filtering can be used 
effectively. 

We discuss a possible stochastic background of gravitational waves 
from the r-modes below. 

Accreting neutron stars (Figure I3.1.4[) are the central objects of 
most of the binary X-ray sources in the Galaxy. Astronomers divide them 
into two distinct groups: the low-mass and high-mass binaries, according 
to the mass of the companion star. In these systems mass is pulled from 
the low or high-mass giant by the tidal forces exerted by its neutron star 
companion. In low-mass X-ray binaries (LMXBs) the accretion lasts long 
enough to spin the neutron star up to the rotation rates of millisecond 
pulsars. Astronomers have therefore supposed for some time that the neu- 
tron stars in LMXBs would have a range of spins, from near zero (young 
systems) to near 500 or 600 Hz (at the end of the accretion phase) . Until 
the launch of the Rossi X-ray Timing Explorer (RXTE), there was no ob- 
servational evidence for the neutron star spins. But in the last two years 
there has been an accumulation of evidence that most, if not all, of these 
stars have angular velocities in a narrow range around 300 Hz.[Tn] It is not 
known yet what mechanism regulates this spin, but a strong candidate is 
the emission of gravitational radiation. 

A novel proposal by Bildsten[20] suggests that the temperature gra- 
dient across a neutron star that is accreting preferentially at its magnetic 
poles should lead to a composition and hence a density gradient in the deep 
crust. Spinning at 300 Hz, such a star could radiate as much as it accretes. 
It would then be a steady source for as long as accretion lasts, which could 
be millions of years. 

In this model the gravitational wave energy flux is proportional to the 
observed X-ray energy flux. The strongest source in this model is Sco X-1, 
which could be detected by GEO600 in a two-year-long narrow-band mode 
if the appropriate matched filtering can be done. LIGO II would have no 
difficulty in detecting this source. 

Older stars may also be lumpy. For known pulsars, this is constrained 
by the rate of spin down: the energy radiated in gravitational waves can not 
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Figure 3.1. Accreting neutron star in a Low-Mass X-Ray binary system 




exceed the total energy loss. In most cases, this limit is rather weak, and 
stars would have to sustain strains in their crust of order 10"'^ or more. It is 
unlikely that crusts could sustain this kind of strain, so the observational 
limits are probably significant overestimates for most pulsars. However, 
millisecond pulsars have much slower spindown rates, and it would be easier 
to account for the strain in their crusts, for example as a remnant Bildsten 
asymmetry. Such stars could in principle be radiating more energy in 
gravitational waves than electromagnetically. 

Observations of individual neutron stars would be rich with informa- 
tion about astrophysics and fundamental nuclear physics. So little is known 
about the physics of these complex objects that the incentive to observe 
their radiation is great. 

However, making such observations presents challenges for data anal- 
ysis, since the motion of the Earth puts a strong phase modulation on 
the signal, which means that even if its rest-frame frequency is constant it 
cannot be found by simple Fourier analysis. More sophisticated pattern- 
matching (matched-filtering) techniques are needed, which track and match 
the signal's phase to within one cycle over the entire period of measure- 
ment. This is not difficult if the source's location and frequency are known, 
but the problem of doing a wide-area search for unknown objects is very 
challenging. [21] Moreover, if the physics of the source is poorly known, such 
as for LMXBs or r-mode spindown, the job of building an accurate family 
of templates is a difficult one. These questions arc the subject of much 
research today, but they will need much more in the future. 
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3.1.5 Random backgrounds 

The big bang was the most violent event of ah, and it may have created 
a significant amount of gravitational radiation. Other events in the early 
universe may also have created radiation, and there may be backgrounds 
from more recent epochs. We have seen earlier, for example, that compact 
binary systems in the Galaxy will merge into a confusion-limited noise 
background in LISA observations below about 1 mHz. 

Let us consider the r-modes as another important example. This pro- 
cess may have occurred in a good fraction of all neutron stars formed since 
the beginning of star formation. The sum of all of their r-mode radiation 
will be a stochastic background, with a spectrum that extends from a lower 
cutoff of about 200 Hz in the rest frame of the emitter to an upper limit 
that depends on the initial angular velocity of stars. If significant star for- 
mation started at, say, a redshift of 5, then this background should extend 
down to about 25 Hz. If 10~^ of the mass of the Galaxy is in neutron 
stars, and each of them radiates 10% of its mass in this radiation, then the 
gravitational wave background should have a density equal to 10"** of the 
mean cosmological density of visible stars. Expressed as a fraction flg^^ of 
the closure density of the Universe, per logarithmic frequency interval, this 
converts to 

^r-modos (25 „ lOOOHz) « 10"^ - 10"^ 

This background would be easily detectable by LIGO II. 

There may also be a cosmological background from either topologi- 
cal defects (e.g. cosmic strings) or from inflation (which amplifies initial 
quantum gravitational fluctuations as it does the scalar ones that lead to 
galaxy formation). Limits from COBE observations suggest that standard 
infiation could not produce a background stronger than V,™^'^^^°'^ ^ 10"^^ 
today. This is too weak for any of the planned detectors to reach, but it 
remains an important long-range goal for the field. But there could also 
be a component of background radiation that depends on what happened 
before infiation: string cosmological models, for example, predict spectra 
growing with frequency. [221 

First-generation interferometers are not likely to detect these back- 
grounds: they may not be able to go below the upper limit set by the 
requirement that gravitational waves should not disturb cosmological nu- 
cleosynthesis, which is ilg^ = 10~^. (This limit does not apply to back- 
grounds generated after nucleosynthesis, like the r-mode background.) Bar 
detectors may do as well or better than the first generation of interfer- 
ometers for a broad-spectrum primordial background: as we have noted 
earlier, their noise levels within their resonance bands are very low. But 
their frequencies are not right for the r-mode background. 

Second-generation interferometers may be able to reach to 10^^^ of 
closure or even lower, by cross-correlation of the output of the two detectors. 
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But they are unlikely to get to the inflation target of 10^^^. LISA may 
be able to go as low as 10^^° (if we have a confldent understanding of the 
instrumental noise) , but it is likely to detect only the confusion background 
of binaries, which is expected to be much stronger than a cosmological 
background in the LISA band. 

3.1.6 The unexpected 

At some level, we are bound to see things we did not expect. LISA, with its 
high signal-to-noise ratios for predicted sources, is particularly well placed 
to do this. Most of the Universe is composed of dark matter whose ex- 
istence we can infer only from its gravitational effects. It would not be 
particularly surprising if a component of this dark matter produced gravi- 
tational radiation in unexpected ways, such as from binaries of small exotic 
compact objects of stellar mass. We will have to wait to see! 



Chapter 4 

Waves and Energy 



Here we discuss wave-like perturbations h^i, of a general background metric 
g^,/. The mathematics is similar to that of linearized theory: /i^^ is a 
tensor with respect to background coordinate transformations (as it was 
for Lorentz transformations in linearized theory) and it undergoes a gauge 
transformation when one makes an infinitesimal coordinate transformation. 
As in linearized theory, we will assume that the amplitude of the waves is 
small. Moreover, the waves must have a wavelength that is short compared 
to the radius of curvature of the background metric. These two assumptions 
allow us to visualize the waves as small ripples running through a curved 
and slowly changing space-time. 



4.1 Variational principle for general relativity 

We start our analysis of small perturbation /i^^ by introducing the standard 
Hilbert variational principle for Einstein's equations. The field equations 
of general relativity can be derived from an action principle using the Ricci 
scalar curvature as the Lagrangian density. The Ricci scalar (second con- 
traction of the Riemann tensor) is an invariant quantity which contains in 
addition to g^i, and its first derivatives also the second derivatives of g^^, 
so our action can be written symbolically as: 

i[gf^,A ^ Ye^ / R{gf,u, g^u.a, g^,u,afi)\f^d'^x (4.1) 

where \/—g is the square-root of the determinant of the metric tensor. As 
usual in variational principles, the metric tensor components are varied 
g^l,l' ^" g^Jl,v + h^i, , and one demands that the resulting change in the action 
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should vanish to first order in any small variation hf^^, of compact support: 

1 



,^ , . -h^,d''x + 0{2) (4.2) 

= " i / ^""'^A^-^^^'-^ + 0(2) (4-3) 

where "0(2)" denotes terms quadratic and higher in hf^^. All the diver- 
gences obtained in the intermediate steps of this calculation integrate to 
zero since /i^i, is of compact support. This variational principle therefore 
yields the vacuum Einstein equations: G^'^ = 0. 

Let us consider how this changes if we include matter. This will help 
us to see how we can treat gravitational waves as a new kind of "matter" 
field on spacetime. 

Suppose we have a matter field, described by a variable $ (which may 
represents a vector, a tensor or a set of tensors). It will have a Lagrangian 
density Lm = Lm (*&, 'I'.q, ■ • ■ l9^J,u) that depends on the field and also on 
the metric. Normally derivatives of the metric tensor do not appear in Lm, 
since by the equivalence principle, matter fields should behave locally as if 
they were in flat space-time, where of course there are no metric derivatives. 
Variations of i,,, with respect to $ will produce the field equation(s) for 
the matter system , but here we are more interested in variations with 
respect to 5^1/, which is how we will find the matter field's contribution to 
the gravitational field equations. The total action has the form: 



/= / {R+l6TTLm)yf^(rx, (4.4) 

whose variation is 



51 = / V V^'^ + / 16^ ^ "^ '" h^.d^'x. (4.5) 
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This variation must yield the full Einstein equations, so we must have the 
following result for the stress-energy tensor of the matter: 

T^-V^^2^^^p^, (4.6) 

leading to 



G>"' = STTTf"" . (4.7) 



This way of deriving the stress-energy tensor of the matter field has deep 
connections to the conservation laws of general relativity, to the way of 
constructing conserved quantities when the metric has symmetries and to 
the so-called pseudo-tensorial definitions of gravitational wave energy (see 
Landau and Lifshitz 1962). [53]. We shall use it in the latter sense. 
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4.2 Variational principles and the energy in gravita- 
tional waves 

Before we introduce the mathematics of gravitational waves, it is impor- 
tant to understand which geometries we are going to examine. We have 
said that these geometries consist of a sfowly and smoothly changing back- 
ground metric which is altered by perturbations of small amplitude and 
high frequency. If L and A are the characteristic lengths over which the 
background and "ripple" metrics change significantly, we assume that the 
ratio ^ will be very much smaller than unity and that \hf^^\ is of the same 
order of smallness as ^ . In this way the total metric remains slowly chang- 
ing on a macroscopic scale, while the high-frequency wave, when averaged 
over several wavelengths, will be the principal source of the curvature of the 
background metric. This is the "short-wave" approximation |24) . Obviously 
this is a direct generalization of the treatment in Lecture 1. 



4.2.1 Gauge transformation and invariance 

Consider an infinitesimal coordinate transformation generated by a vector 
field ^", 

x^-^x^+C- (4.8) 

In the new coordinate system, neglecting quadratic and higher terms in 
/i"^, it is not hard to show that the general gauge transformation of the 
metric is 

where a semicolon denotes the covariant derivative. We assume that the 
derivatives of the coordinate displacement field are of the same order as 
the metric perturbation: |^"'^| ~ Ift."'^!. 

Isaacson 24 showed that the gauge transformation of the Ricci and 
Riemann curvature tensors has the property 

Rl'J-Rl'J^ (^)' (4.10) 

where Rl,J and -R^j^^, are the first order Ricci and Riemann tensors (in 
powers of perturbation h^^) and an over-bar denotes their values after the 
gauge transformation. In our high-frequency limit, therefore, these tensors 
are gauge-invariant to linear order, just as in linearized theory. 
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4.2.2 Graviational wave action 

Let us suppose that we are in vacuum so we iiave only tlie metric, no matter 
fields, but we work in the high-frequency approximation. The full metric is 
g^^ (smooth background metric) +h^i, (high-frequency perturbation) . Our 
purpose is to show that the wave field can be treated as a "matter" field, 
with a Lagrangian and its own stress-energy tensor. To do this we have to 
expand the action out to second order in the metric perturbation. 



I[gfifi + h^u] = / R{gij.u + h^^, gi_,y^a + /ipi/.a, ■ • ■)\l -g[9tiu + h^i,]d'^x 

' r ^ ( T? I — ^ 
R{3ii.v, ■ ■ ■)y^d*x + / — hf.^d'^x 



1 ffd^RV^) ^2^-1:^^^/, h 



2 y V dg^^^dgap dg^j^^dgaij,^ 



dgf_Lu,rdgai3,'f ' ' dgi_,^dga/3,fT 

+ 0(3). 

The first term is the action for the background metric 5^,^. The second 
term vanishes (see Equation (J4.2I) '). since we assume that the background 
metric is a solution of Einstein vacuum equation itself, at least to lowest 
order. 

If we compare the above equation with Equation (14. 4L we can see 
that the third term, complicated as it seems, is an effective "matter" La- 
grangian for the gravitational field. Indeed, if one varies it with respect 
to hfii, holding g^^, fixed (as we would do for a physical matter field on 
the background), then the complicated coefficients are fixed and one can 
straightforwardly show that one gets exactly the linear perturbation of the 
Einstein tensor itself. Its vanishing is the equation for the gravitational 
wave perturbation /i^^. In this way we have shown that, for a small am- 
plitude perturbation, the gravitational wave can be treated as a "matter" 
field with its own Lagrangian and field equations. 

Given this Lagrangian, we should be able to calculate the effective 
stress-energy tensor of the wave field by taking the variations of the effective 
Lagrangian with respect to gf^i,, holding the "matter" field /i^i, fixed: 

j^(GW)a0 / — ^=2 ^'^ [gMi^ifeA'i^]\/~g (4 11) 

dgai3 
with 

327r V og^^dgafs ^ "^ ogf,uOgap,y 
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This quantity is quadratic in the wave amphtude ft,^^. It could be simplified 
further by integrations by parts, such as by taking a derivative off hap^-yr- 
This would change the coefficients of the other terms. We will not need to 
worry about finding the "best" form for the expression (4.12), as we now 
show. 

As in linearized theory, so also in the general case, the quantity /i^^ 
behaves like a tensor with respect to background coordinate transforma- 
tions, and so does T/;^ . But it is not gauge-invariant and so it is not a 
physical observable. However, since the integral of the action is indepen- 
dent of coordinate transformations that have compact support, so too is the 
integral of the effective stress-energy tensor. In practical terms, this makes 
it possible to localize the energy of a wave to within a region of about one 
wavelength in size where the background curvature does not change signif- 
icantly. In fact, if we restrict our gauge transformations to have a length 
scale of a wavelength, and if we average (integrate) the stress-energy ten- 
sor of the waves over such a region, then any gauge changes will be small 
surface terms. 

By evaluating the effective stress-energy tensor on a smooth back- 
ground metric in a Lorcntz gauge, and performing the averaging (denoted 
by symbol (• • • ) ), one arrives at the Isaacson tensor: 



^.(GW) 



^ ( V;^/^'' ;^) • (4-13) 
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This is a convenient and compact form for the gravitational stress- 
energy tensor. It localizes energy in short-wavelength gravitational waves 
to regions of the order of a wavelength. It is interesting to remind ourselves 
that our only experimental evidence of gravitational waves today is the 
observation of the effect on a binary orbit of the loss of energy to the 
gravitational waves emitted by the system. So this energy formula, or 
equivalent ones, is central to our understanding of gravitational waves. 

4.3 Practical applications of the Isaacson energy 

If we are far from a source of graviational waves, we can treat the waves 
by linearized theory. Then if we adopt the TT gauge and specialize the 
stress-energy tensor of the radiation to a flat background, we get 

Since there are only two components, a wave travelling with frequencies 
/ (wave number k = 27r/) and with a typical amplitude h in both po- 
larizations carries an energy Fg^ equal to (see Ex.. 6 at the end of this 
lecture) 

Fg^ = jfh\ (4.15) 
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Putting in the factors of c and G and scaling to reasonable values gives 
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(4.16) 



which is a very large energy flux even for this weak a wave. It is twice the 
energy flux of a full moon! Integrating over a sphere of radius r, assuming 
a total duration of the event r, and solving for h, again with appropriate 
normalizations, gives 
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(4.17) 



This is the formula for the "burst energy" , normalized to numbers appro- 
priate to a gravitational collapse occurring in the Virgo cluster. It explains 
why physicists and astronomers regard the 10 ~^^ threshold as so impor- 
tant. But this formula could also be applied to a binary system radiating 
away its orbital gravitational binding energy over a long period of time r, 
for example. 

4.3.1 Curvature produced by waves 

W have assumed that the background metric satisfied the vacuum Einstein 
equations to linear order, but now it is possible to view the full action 
principles as a principle for the background with a wave field h^i, on it, and 
to let the wave energy affect the background curvature. [23 This means that 
the background will actually solve, in a self-consistent way, the equation 



Gap [g^iu] 






[9^l^l^ 



(4.18) 



This does not contradict the vanishing of the first variation of the action, 
which we needed to use above, because now we have an Einstein tensor 
that is of quadratic order in /i^^, contributing a term of cubic order to the 
first-variation of the action, which is of the same order as other terms we 
have neglected. 



4.3.2 Cosmological background of radiation 

This self-consistent picture allows us to talk about, for example, a cosmo- 
logical gravitational wave background that contributes to the curvature of 
the Universe. Since the energy density is the same as the flux (when c = 1), 
we have 

Pgw = ^fh', (4.19) 

but now we must interpret h in a statistical way. This will be treated in 
the contribution by Babusci et al., but basically it is done by replacing h^ 
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by a statistical mean square amplitude per unit frequency (Fourier trans- 
form power), so that the energy density per unit frequency is proportional 

to f^ h . It is then conventional to talk about the energy density per 
unit logarithm of frequency, which means multiplying by /. The result, 
after being careful about averaging over all directions of the waves and all 
independent polarization components, is 






^47r'f\h{f)\\ (4.20) 



Finally, what is of the most interest is the energy density as a fraction of 
the closure or critical cosmological density, given by the Hubble constant 
Hq as Qc = 3Hq/8tt. The resulting ratio is the symbol r2gw(/) that we met 
in the previous Lecture: 
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^sM) = ^f\hif)\. (4.21) 



4.3.3 Other approaches 

We finish this Lecture by observing that there is no unique approach to 
defining energy for gravitational radiation or indeed for any solution of 
Einstein's equations. Historically this has been one of the most difficult 
areas for physicists to come to grips with. In the textbooks you will find 
discussions of pseudotensors, of energy measured at null infinity and at 
spacelike infinity, of Noether theorems and formulas for energy, and so on. 
None of these are worse than we have presented here, and in fact all of 
them are now known to be consistent with one another, if one does not ask 
them to do too much. In particular, if one wants only to localize the energy 
of a gravitational wave to a region of the size of a wavelength, and if the 
waves have short wavelength compared to the background curvature scale, 
then pseudotensors will give the same energy as the one we have defined 
here. Similarly, if one takes the energy flux defined here and evaluates it 
at null infinity, one gets the so-called Bondi flux, which was derived by H. 
Bondi in one of the pioneering steps in the understanding of gravitational 
radiation. Many of these issues are discussed in the Schutz-Sorkin paper 
referred to earlier. [23] 



4.4 Exercises 

5 In the notes above we give the general gauge transformation 

i^^u ' l^fiu >,^:y ^^;/^' 
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Use the formula for the derivation of Einstein's equations from an 
action principle, 

1 fSiRV^j)^ . 



with 

Ur Y i/l 



but insert a pure gauge h^^. Argue that since this is merely a coor- 
dinate transformation, the action should be invariant. Integrate the 
variation of the action to prove the contacted Bianchi identity 

Gf\^ ^ 0. 

This shows that the divergence-free property of G^'^ is closely related 
to the coordinate invariance of Einstein's theory. 

6 Suppose a plane wave, travelling in the z-direction in linearized theory, 
has both polarization components /i+ and hx • Show that its energy flux 
in the z-direction, r(GW)Oz^ ^^ 

where the angle brackets denote an average over one period of the wave. 



Chapter 5 

Mass- and Current-Quadrupole 
Radiation 



In this Lecture we focus on the wave ampUtude itself, and how it and the 
polarization depend on the motions in the source. Consider an isolated 
source with a stress-energy tensor T"^. As in Lecture 1, the Einstein 
equation is 

^' ■ V^) r" = -16^r"'5 (5.1) 



9t2 

{h — h"^ — ^■q°'^h and h « = 0). Its general solution is the following 
retarded integral for the field at a position x* and a time t in terms of the 
source at a position y* and the retarded time t — R: 

r^ {x\t)=4j |t"^ {t - R, y^) £y, (5.2) 

where we define 

i?2 = (x^ - y^) (x, - y,) . (5.3) 

5.1 Expansion for the far field of a slow-motion source 

Let us suppose that the origin of coordinates is in or near the source, and 
the field point a;* is far away. Then wc define r^ = x'^Xi and wc have 
3> y^yi- We can therefore expand the term R in the dominator in terms 



r 



of y*. The lowest order is r, and all higher-order terms are smaller than 
this by powers of r^^. Therefore, they contribute terms to the field that 
fall off faster than r~^, and they are negligible in the far zone. So we can 
simply replace i? by r in the dominator, and take it out of the integral. 

The R inside the time-argument of the source term is not so simple. 
If we suppose that T"^ does not change very fast we can substitute t — R 
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by t — r (the retarded tinie to the origin of coordinates) and expand 

t- R = t-r + n^yi + { -] , with n* = — , nVj = 1. (5.4) 

\r J r 

The two conditions r ^ y^yi and the slow-motion source, can be ex- 
pressed quantitatively as: 

r > A 

i?<A 

where A is the reduced wave length A — A/27r and R is the size of source. 
The terms of order r^^ are negligible for the same reason as above, 
but the first term in this expansion must be taken into account. It depends 
on the direction to the field point, given by the unit vector n*. We use this 
by making a Taylor expansion in time on the time-argument of the source. 
The combined effect of these approximation is 



-a/3 4 

h ~ — 



r _ 
1 



T"^ (f, y^) + T"% (f, f) n^y, + \t"^^,, (f, f) n=n^y,yu 

£y. (5.5) 



We will need all the terms of this Taylor expansion out to this order. 

The integrals in expression (5.5) contain moments of the components 
of the stress-energy. It is useful to give these names. Use M for moments 
of the density T*"^, P for moments of the momentum T°*, and S for the 
moments of the stress T*^. Here is our notation: 



M (f) = J T°" (f, f) £y, M' (f) = j T°" (f, f) y^d\ 
^,k (^,) ^ f yoo (^,^ y^^ y^y^d'y, M=^' (f) = f T^^ {f , f) y^y'y^d\ 



P' it') = / T°' (i', f) d\ P'-'" it') = / T°' {t', y^) y'd^y, 



pm^t')^ T^'{t\f)y^y''d'y, 



glm ^^r^ ^ / rj^lm ^, ^ y,^ ^3^^ glmj ^^r^ ^ / rj.lm ^, ^ y^^ yj^3y_ 



These are the moments we will need. 

Among these moments there are some identities that follow from the 

,/3 



conservation law in linearized theory, T" „ — 0, which we use to replace 
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time derivatives of components of T by divergences of other components 
and then integrate by parts. The identities we will need are 



pj ^ 0, P''' = S'^■^ P'''^ = S'''^ + S'^^. (5.7) 

These can be applied recursively to show, for example, two further very 
useful relations 

where the round brackets on indices indicate full symmetrization. 
Using this relations and notations it is not hard to show that 

Tf\t,x') = -M + -P^rij + -S'''{t')njnk + -S^''\t')njnkni + . . (5.9) 

TfUt,x') = -P^ + -S^''(t')nk + -S^''\t')nkni + ... (5.10) 

r r r 

h'\t,x') = -S^'^it') + -S'''\t')ni + .... (5.11) 

r r 

In these three formulas there are different orders of time-derivatives, but 
in fact they are evaluated to the same final order in the slow-motion ap- 
proximation. One can see that from the gauge condition h « = 0, which 
relates time-derivatives of some components to space-derivatives of others. 
In these expressions, one must remember that the moments are eval- 
uated at the retarded time t' — t — r (except for those moments that are 
constant in time), and they are multiplied by components of the unit vector 
to the field point n^ — x^ j r. 



5.2 Application of TT gauge to the mass quadrupole 
field. 

In the expression for the amplitude that we derived so far, the final terms 
are those that represent the current-quadrupole and mass-octupole radi- 
ation. The terms before them represent the static parts of the field and 
the mass quadrupole radiation. In this section we treat just these terms, 
placing them into TT gauge. This will be simpler than treating it all at 
once, and the procedure for the next terms will be a straightforward gen- 
eralization. 
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5.2.1 The TT gauge transformations 

We are already in Lorentz gauge, and this can be checked by taking deriva- 
tives of the expressions for the field that we have derived above. But we 
are manifestly not in TT gauge. Making a gauge transformation consists 
of choosing a vector field ^" and modifying the metric by 

haf! ^ haj3 — S,a,l3 — ^I3,a- (5-12) 

—a/3 

The corresponding expression for the potential h is 

r^^r%r'^+r'^-r?"''?%. (5.13) 



A* 



For the different components this implies changes 



<5r =e°'°+e^, (5.14) 

S7f^ =^^-3 +^3'" (5.15) 



Sh^'" = e^'*^ + e'^^^' - 5^''^,^ (5.16) 

where S^'^ is the Kronecker delta (unit matrix). In practice, when taking 
derivatives, the algebra is vastly simplified by the fact that we are keeping 
only r~^ terms in the potentials. This means that spatial derivatives do 
not act on r^^ but only on i' := t — r. It follows that dt' /dx^ = —Uj, and 
dh{t')/dx3 = -h{t')nj. 

It is not hard to show that the following vector field puts the metric 
into TT gauge to the order we are working: 

^° = ip\ + -P'''n,nk + -S\kn'' + -S'^^'n.n.nk, (5.17) 

r r r r 

e = -M' + -P'^Ui - -P\n' - -P^^n.nku' + -S'^'^n^nu 

- -S\kn''n' - -S'^^njUinunK (5.18) 

5.2.2 Quadrupole field in TT gauge 

The result of applying this gauge transformation to the original amplitudes 
is 

7^™° = ^, (5.19) 

r 

7^^™ = 0, (5.20) 



t-TTjj _ 4 



V'^l.^' Sik + \ r^ {Skin'^n' - ^\) 



(5.21) 
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Remember that here wc are not including S'^'^', because it is a third order 
effect. 

The notation 1. represents the projection operator perpendicular to 
the direction n* to the field point. 

j_jk^ pk _ ^3^k_ (5 22) 

It can be verified that this tensor is transverse to the direction n* and is a 
projection, in the sense that it projects to itself 

L^^ Uk = 0, _L^''=_Lfc'=_L^' . (5.23) 

The spherical component of the field is not totally eliminated in this 
gauge transformation: the time-time component of the metric must contain 
the constant Newtonian field of the source. (In fact we have succeeded in 
eliminating the dipole, or momentum part of the field, which is also part 
of the non-wave solution. Our gauge transformation has incorporated a 
Lorentz transformation that has put us into the rest frame of the source.) 
The time-dependent part of the field is now purely spatial, transverse (be- 
cause everything is multiplied by _L), and traceless (as can be verified by 
explicit calculation). 

The expression for the spatial part of the field actually does not depend 
on the trace of Sjk, as can be seen by constructing the trace-free part of 
the tensor, defined as: 

SJk ^ g]k _ hjkgi^ (5 24) 

In fact, it is more conventional to use the mass moment here instead of the 
stress, so we also define 

M3^ ^M^^ --5^^Ml S^^^-——-. (5.25) 

In terms of M the far field is 

h^^'' = ^ U'^^^'Mki +i ±''Mki n'n''^ . (5.26) 

This is the usual formula for the mass-quadrupole field. In textbooks the 
notation is somewhat different than we have adopted here. In particular, 
our tensor M is what is called I in Misner, et al, (1973) and Schutz (1985). 
It is the basis of most gravitational wave source estimates. We have derived 
it only in the context of linearized theory, but remarkably its form is identi- 
cal if we go to the post- Newtonian approximation, where the gravitational 
waves are a perturbation of the Newtonian spacetime rather than of flat 
spacetime. 
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Given this powerful formula, it is important to try to interpret it and 
understand it as fully as possible. One obvious conclusion is that the dom- 
inant source of radiation, at least in the slow- motion limit, is the second 
time-derivative of the second moment of the mass density T^'^ (the mass- 
quadrupole moment). This is a very important difference between gravita- 
tional waves and electromagnetism, in which the most important source is 
the electric-dipole. In our case the mass-dipole term is not able to radiate 
because it is constant, reflecting conservation of the linear momentum of 
the source. In electromagnetism, however, if the dipole term is absent for 
some reason (all charges positive, for example) then the quadrupole term 
dominates and it looks very similar to Equation (j5.26p . 

5.2.3 Radiation patterns related to the motion of sources 

The projection operators in Equation (j5.26p show that the radiative field 
is transverse, as we expect. But the form of Equation (|5.26p hides two 
equally important messages, 

• the only motions that produce the radiation are the ones transverse 
to the line of sight; and 

• the induced motions in a detector mirror the motions of the source 
projected onto the plane of the sky. 

To see why these are true, we define the transverse traceless quadrupole 
tensor 

Mf/' =± \ ± 'jMki - ^ ±„a'=' Mki. (5.27) 

(Notice that some of our definitions of tracelessness involve subtracting 1/3 
of the trace, as in Equation (j5.24|) . and sometimes 1/2 of the trace, as here 
in Equation (J5.27I) . The appropriate factor is determined by the effective 
dimensionality (rank) of the tensor. Although we have 3 spatial dimen- 
sions, the projection tensor _L projects the mass quadrupole tensor onto a 
two-dimensional plane, where the trace involves only two components, not 
three.) 

Now, if in Equation (j5.26p we replace Mij by its definition in terms of 
Mij, and then collect terms appropriately, it is not hard to show that the 
equation simplifies to its most natural form: 

-h'^^^J = -M™'. (5.28) 

r 

This could of course have been derived directly by applying the TT oper- 
ation to Equation (J5.9I) to 15.111 Since this equation involves only the TT- 
part of M, our first assertion above is proved. According to this equation, 
in order to calculate the quadrupole radiation that a particular observer 
will receive, one need only compute the mass-quadrupole tensor's second 
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time-derivative, project it onto the plane of the sky as seen by the observer 
looking toward the source, take away its trace, and rescale it by a factor 
of 2/r. In particular, the TT-tensor that describes the action of the wave 
(as in the polarization diagram in Figure II. ip is a copy of the TT-tensor of 
the mass distribution. This proves our second assertion above. 

Look again at Figure 11.11 Imagine a detector consisting of two free 
masses whose separation is being monitored. If the wave causes them to 
oscillate relative to one another along the x-axis (the ® polarization), this 
means that the source motion contained a component that did the same 
thing. If the source is a binary, then the binary orbit projected onto the 
sky must involve motion of the stars back and forth along either the x- or 
the y-axis. 

It is possible from this to understand many aspects of quadrupole 
radiation in a simple way. Consider a binary star system with a circular 
orbit. Seen by a distant observer in the orbital plane, the projected source 
motion is linear, back and forth. The received polarization will be linear, 
the polarization ellipse aligned with the orbit. Seen by a distant observer 
along the axis of the orbit of the binary, the projected motion is circular, 
which is a superposition of two linear motions separated in phase by 90°. 
The received radiation will also have circular polarization. Because both 
linear polarizations are present, the amplitude of the wave emitted up the 
axis is twice that emitted in the plane. In this way we can completely 
determine the radiation pattern of a binary system. 

Notice that, when viewed at an arbitrary angle to the axis, the radia- 
tion will be elliptically polarized, and the degree of ellipticity will directly 
measure the inclination of the orbital plane to the line of sight. This is a 
very special kind of information, which one cannot normally obtain from 
electromagnetic observations of binaries. It illustrates the complementarity 
of the two kinds of observing. 



5.3 Application of TT gauge to the current quadrupole 
field 



Now we turn to the problem of placing next-order terms of the wave field, 
the current-quadrupole and mass-octupole, into TT-gauge. Our interest 
here is to understand current-quadrupole radiation in the same physical 
way as we have just done for mass-quadrupole radiation. So we shall put 
the field into TT gauge and then see how to separate the current-quadrupole 
part from the mass-octupole, which we will discard from the present dis- 
cussion. 
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5.3.1 The field at third order in slow-motion 

The next order terms in the non-TT metric bear a simple relationship to 
the mass quadrupole terms (see equations 15. 9115. iTj) . In each of the metric 
components, just replace S^'^ by S^^^ni to go from one order to the next. 

This means that we can just skip to the end of the application of the 
gauge transformations in equations 15.171 and 15.181 and write the next order 
of the final field, only using S again, not M: 

h^^'' = J ±''±^' ^ifcrnn" + i ±'' (^fci^n^n'n™ - ^\;n') , (5.29) 

or more compactly 

h™' = J (±"'±''Skim n"' + i r^Skim n^n'^uA . (5.30) 

The tilde above S represents a trace-free operation on the first two indices, 

These are the indices that come from the indices of T^'^, so the tensor is 
symmetric on these. By analogy with the quadrupole calculation, we can 
also define the TT part of Sijk by doing the TT projection on the first two 
indices, 

^ijm ^-'- i -'- j^klm ^ -^ij-^ i'klm- (O.olj 

The TT projection of the equation for the metric is 

r 

5.3.2 Separating the current-quadrupole from the mass-octupole 

The last equation is compact, but it does not have the ready interpreta- 
tion that we have at quadrupole order. This is because the moment of 
the stress, Sijk, does not have such a clear physical interpretation. We see 
from Equation (j5.6p to 15.81 that Sijk is a complicated mixture of moments 
of momentum and density. To gain more physical insight into radiation at 
this order, we need to separate these different contributions. It is straight- 
forward algebra to see that the following identity follows from the earlier 
ones: 

S'^f" = - M '^^ + - P [^■''1' + - P ['''1^ (5.33) 

6 3 3 

where square brackets around indices mean antisymmetrization: 

j^[ik] ._ £ I'^ife __ j^ki\ 
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This is a complete separation of the mass terms (in M) from the momentum 
terms (in P) because the only identities relating the momentum moments 
to the mass moments involve the symmetric part of P*^*^ on its first two 
indices, and this is absent from Equation (j5.33p . 

The first term in Equation (|5.33p is the third moment of the density, 
and this is the source of the mass octupole field. It produces radiation 
through the third time-derivative. Since we are in a slow-motion approxi- 
mation, this is smaller than the mass quadrupole radiation by typically a 
factor of v/c. Unless there were some very special symmetry conditions, 
one would not expect the mass octupole to be anything more than a small 
correction to the mass quadrupole. For this reason we will not treat it here. 

The second and third terms in Equation (j5.33p involve the second 
moment of the momentum, and together they are the source of the cur- 
rent quadrupole field. It involves two time-derivatives, just as the mass- 
quadrupole does, but these are time-derivatives of the mometum moment, 
not the mass moment, so these terms produce a field that is also v/c smaller 
than the typical mass quadrupole field. However, it requires less of an ac- 
cident for the mass quadrupole to be absent and the current quadrupole 
present. It just requires motions that leave the density unchanged to lowest 
order. This happens in the r-modes. Therefore, the current-quadrupole de- 
serves more attention, and we will work exclusively with these terms from 
now on. 

The terms in Equation (|5.33p that we need are the ones involving 

P *•'*. These are atni-symmetrized on the first two indices, which involves 
effectively a vector product between the momentum density (first index) 
and one of the moment indices. This is essentially the angular momentum 
density. To make the angular momentum explicit and to simplify the ex- 
pression, we introduce the angular momentum and the first moment of the 
angular momentum density 

J' := e'^'^Pjk, (5.34) 

r' := e'^'^Pjfc', (5.35) 

where e^^'^ is the fully antisymmetric (Levi-Civita) symbol in 3 dimensions. 
It follows from this that 

2 
These terms enter the TT projection of the field Equation ()5.32p with the 
last index of S always contracted with the direction n^ to the observer from 
the source. According to Equation (j5.33p . this contraction always occurs 
on one of the antisytmmetrized indices, or if we use the form in the previous 
equation then we will always have a contraction of rf with e^^^ . This is a 
simple object, which we call 

^e'^ := n,e'^^. (5.36) 
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This is just the two-dimensional Levi-Civita object in the plane perpendic- 
ular to rf, which is the plane of the sky as seen by the observer. These 
quantities will be used in the currcnt-quadrupole field, which contains pro- 
jections on all the indices. Therefore, the only components of J^'' that 
enter are those projected onto the sky, and so it will simplify formulas to 
define the sky-projected moment of the angular momentum ±J 

iJ*-'":=±^±J„j'™. (5.37) 

Using this assembled notation, the current-quadrupole field is 

^™' = ^ (i^'" ±Jk' + ±6'" ± J k'+ 1-'' ±e'" ±Jkm). (5.38) 

This is similar in form and complexity to the mass-quadrupole field 
expression. The interpretation of the contributions is direct. Only compo- 
nents of the angular momentum in the plane of the sky contribute to the 
field. Similarly only moments of this angular momentum transverse to the 
line of sight contribute. If one wants, say, the xx component of the field, 
then the ±e factor tells us it is determined by the y-componcnt of momen- 
tum, ie the component perpendicular to the x-direction in the sky. In fact, 
it is much simpler just to write out the actual components, assuming that 
the wave travels toward the observer along the z-dircction. Then we have 

/^TTxx ^ 4 ^ ,^^ -^ ^ (5.39) 

6r 

h'^^y^^Qyy-j^^), (5.40) 

6r 

and the remaining components can be found from the usual symmetries of 
the TT-metric. I have dropped the prefix _l on J because in this coordinate 
system the given components are already transverse. 

The simplicity of these expressions is striking. There are two basic 
cases where one gets current-quadrupole radiation. 

• If there is an oscillating angular momentum distribution with a dipole 
moment along the angular momentum axis, as projected onto the sky, 

then in an appropriate coordinate system J ^^ will be nonzero and we 
will have ® radiation. To have a non- vanishing dipole moment, the 
angular momentum density could, for example, be symmetrical under 
reflection through the origin along its axis, so that it points in opposite 
directions on opposite sides. 

• If there is an oscillating angular momentum distribution with a dipole 
moment along an axis perpendicular to the angular momentum axis, 
as projected onto the sky, then in an appropriate coordiate system 

J ^y will be nonzero and we will have ® radiation. 
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5.3.3 A model system radiating current-quadrupole radiation 

To see that the first of these two leads to physicahy sensible results, let 
us consider a simple model system that actually bears a close resemblance 
to the r-mode system. Imagine, as in the left panel of Figure 15.11 two 
wheels connected by an axis, and the wheels are sprung on the axis in such 
a way that if a wheel is turned by some angle and then released, it will 
oscillate back and forth about the axis. Set the two wheels into oscillation 
with opposite phases, so that when one wheel rotates clockwise, the other 
rotates counterclockwise, as seen along the axis. 



Figure 5.1. A simple current-quadrupole radiator. The left panel shows how 
the two wheels are connected with blade springs to a central axis. The wheels 
turn in opposite directions, each oscillating back and forth about its rest position. 
The right panel shows the side view of the system, and the arrows indicate the 
motion of the near side of the wheels at the time of viewing. The + signs indicate 
where the momentum of the mass of the wheel is toward the viewer and the - 
signs indicate where it is away from the viewer. 
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Then when viewed along the axis, the angular momentum has no com- 
ponent transverse to the line of sight, so there is no radiation along the axis. 
This is sensible, because when projected onto the plane of the sky the two 
wheels are performing exactly opposite motions, so the net effect is that 
there is zero projected momentum density. 

When viewed from a direction perpendicular to the axis, with the axis 
along the x-direction, then the angular momentum is transverse, and it has 
opposite direction for the two wheels. There is therefore an a;-moment of 
the x-component of angular momentum, and the radiation field will have 
the igD orientation. 

To see that this has a physically sensible interpretation, look back 
again at the polarization diagram. Figure II. 1[ and look at the bottom row 
of figures illustrating the (g) polarization. See what the particles on the 
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a;-axis are doing. They are moving up and down in the y-direction. What 
motions in the source could be producing this? 

At first one might guess that it is the up-and-down motion of the mass 
in the wheels as they oscillate, because in fact the near side of each wheel 
does exactly what the test particles at the observer are doing. But this 
cannot be the explanation, because the far side of each wheel is doing the 
opposite, and when they both project onto the sky they cancel. What in 
fact gives the effect is that at the top of the wheel the momentum density is 
first positive (towards the observer) and then negative, while at the bottom 
of the wheel it is first negative and then positive. On the other wheel, the 
signs are reversed. 

Current quadrupole radiation is produced, at least in simple situations 
like the one we illustrate here, by the (second time- derivative of) the com- 
ponent of source momentum along the line of sight. If this is positive in the 
sense that it is towards the observer, then the momentum density acts as a 
positive gravitational "charge". If negative, then it is a negative "charge". 
The wheels have an array of positive and negative spots that oscillates with 
time, and the test particles in the polarization diagram are drawn toward 
the positive ones and pushed away from the negative ones. Interestingly, 
in electromagnetism, magnetic dipole and magnetic quadrupole radiation 
are also generated by the component of the electric current along the line 
of sight. 

This is a rather simple physical interpretation of some rather more 
complex equations. It is possible to re- write Equation (j5.38p to show ex- 
plicitly the contribution of the line-of-sight momentum, but the expressions 
become even more complicated. Instead of dwelling on this, I will turn to 
the question of calculating the total energy radiated by the source. 

5.4 Energy radiated in gravitational Waves 

We have calculated the energy flux in Equation (J4.14L and we now have 
the TT wave amplitudes. We need only integrate the flux over a distant 
sphere to get the total luminosity. We do this for the mass and current 
quadrupoles in separate sections. 

5.4.1 Mass quadrupole radiation 

The mass quadrupole radiation field in Equation (J5.26I) must be put into 
the energy flux formula, and the dependence on the direction n* can then be 
integrated over a sphere. It is not a difficult calculation, but it does require 
some angular integrals over over multiple products of the vector n* , which 
depends on the angular direction on the sphere. By symmetry, integrals 
of odd numbers of factors of n* vanish. For even numbers of factors, the 
result is essentially determined by the requirement that after integration 
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the result must be fully symmetric under interchange of any two indices 
and it cannot have any special directions (so it must depend only on the 
Kronecker delta S^j). The identities we need are 

^'n^dn^—S'^, (5.41) 

n'n^n^n^dn = — (5'^ 5^^ + 5'''5^^ + S'^S^'') . (5.42) 

15 

Using these, one gets the following simple formula for the total luminosity 
of mass-quadrupole radiation 



L™"- = i /uJ'^M.k 



(5.43) 



Here we still preserve the angle brackets of Equation (|4.14p . because this 
formula only makes sense in general if we average in time over one cycle of 
the radiation. 

5.4.2 Current quadrupole radiation 

The energy radiated in the current quadrupole is nearly as simple to obtain 
as the mass quadrupole formula. The extra factor of n* in the radiation 
field makes the angular integrals longer, and requires two further identities: 

n'n^n^n^nPn'idVt ^ 5(^1 6^^ 5P1) ^ (5.44) 

^ijk^i'j'k' _ gii' gjj' gkk' _|_ gij' gjk' ^ki' ^ gik' ^ji' gkj' 

- 5'' 5^^' 5^^' - S'^' S^'' d''''' - S'''' 6^^' S'"' , (5.45) 

where the round brackets indicate full symmetrization on all indices. The 
expression is simplest if we define 

jik . ^ilm r> k , klm t3 7 

where 

pkij ._ pkij _ _xij pkl 

3 

The result of the integration of the flux formula over a distant sphere 
is[25][T8|, in our notation, 



j^current ^ _ ^ jjkj^^ ^ (5 ^g) 
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5.5 Radiation in the Newtonian limit 

The calculation so far has been within the assumptions of linearized the- 
ory. Real sources are likely to have significant self-gravity. This means, in 
particular, that there will be a significant component of the source energy 
in gravitational potential energy, and this must be taken into account. 
In fact a more realistic equation than Equation (j5.ip would be 

nh°'^ = -IGtt {T°'^ + t"!^^) (5.47) 

where t"^ is the stress-energy pseudotensor of gravitational waves. This 
is hard to work with: Equation (I5.47P is an implicit equation because i"^ 
depends on h"^ . 

Fortunately, the formulas that we have derived are more robust than 
they seem. It turns out that the leading order radiation field from a New- 
tonian source has the same formula as in linearized theory. By leading 
order we mean the dominant radiation. If there is mass-quadrupole radia- 
tion, then the mass-octupole radiation from a Newtonian source will not be 
given by the formulas of the linearized theory. On the other hand, current 
quadrupole and mass quadrupole radiation can co-exist, because they have 
different symmetries, so the work we have done here is generally applicable. 

More details on how one calculates radiation to higher order in the 
Newtonian limit will be given in the lecture by Blanchet. This is particu- 
larly important for computing the radiation to be expected from coalescing 
binary systems, whose orbits become highly relativistic just before coales- 
cence and which are therefore not well described by linearized theory. 



Chapter 6 

Source calculations 



Now that we have the formulas for the radiation from a system, we can use 
them for some simple examples. 

6.1 Radiation from a binary system 

The most numerous sources of gravitational waves are binary stars systems. 
In just half an orbital period, the non-spherical part of the mass distribution 
returns to its original configuration, so the angular frequency of the emitted 
gravitational waves is twice the orbital angular frequency. 

We shall calculate here the mass quadrupole moment for two stars 
of masses mi and ?7i2, orbiting in the x-y plane in a circular orbit with 
angular velocity il, governed by Newtonian dynamics. We take their total 
separation to be R, which means that the orbital radius of mass mi is 
m,2R/(mi + 7712) while that of mass m2 is miR/{mi + 7712). We place the 
origin of coordinates at the center of mass of the system. Then for example 
the xx-component of M^^ is 

'r7i2-Rcos(r2i)\ f miRcos(nt)\ 

mi I + 7712 ' 



mi + 1712 ) V "^1 + ^^2 / 

= ^i?2^os2(m), (6.1) 

where /z :— mim2/{mi+m,2) is the reduced mass. By using a trigonometric 
identity and throwing away the part that does not depend on time (since 
we will use only time-derivatives of this expression) we have 

M^^ = -fiR'^cos{2nt). (6.2) 

By the same methods, the other non-zero components are 

Myy = — ^iR^ cos{2nt), M^y ^ -^iR^ sin(2m). 
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This shows that the radiation will come out at twice the orbital frequency. 

In this case the trace-free moment M*-' differs from M *■' only by a 
constant, so we can use these values for M*^ to calculate the field and 
luminosity. 

As an example of calculating the field, let us compute /iTTxx g^g gg^j^ 
by an observer at a distance r from the system along the y-axis, i.e. lying 
in the plane of the orbit. We first need the TT-part of the mass quadrupole 
moment, from Equation (|5.27p : 

But since Af "^^ = 0, this is just M'==^/2. Then from Equation ((5?28l) we find 

jjTxx ^ _2^(i?f])2 cos [2n{t - r)] . (6.3) 

r 



Similarly, the result for the luminosity is 

32 



Lgu, - — M'i?^^'- (6.4) 



The various factors in thse two equations are not independent, be- 
cause the angular velocity is determined by the masses and separations of 
the stars. When observing such a system, we can not usually measure R 
directly, but we can infer n from the observed gravitational wave frequency, 
and we may often be able to make a guess at the masses. (Actually, we 
will see below that we can actually measure the important quantity about 
the masses.) So we eliminate R using the Newtonian orbit equation 

R^ - ^^^. (6.5) 

If in addition we use the gravitational wave frequency ilg^ — 251, we get 

-^TTxx ^ _ 21/3^^:1_1Vl COS pg^{t - r)] , (6.6) 

Lg^ = ^^{Mn,^)'^ , (6.7) 

where we have introduced the symbol for the chirp mass of the binary 
system: 

Notice that both the field and the luminosity depend only on A^, not on 
the individual masses in any other combination. 

The power represented by L^u, must be supplied by the orbital energy, 
E = —niim2/2R. By eliminating R as before we find the equation 

^^ 25/3 "ff"^ ■ 
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This is remarkable because it too involves only the chirp mass Ai. By 
setting the rate of change of E equal to the (negative of the) luminosity, we 
find an equation for the rate of change of the gravitational wave frequency 

n,^ = ^-^^M'/'nli/\ (6.8) 

As we mentioned in Lecture 3, since the frequency increases, the signal is 
said to "chirp" . 

These results show that the chirp mass is the only mass associated 
with the binary that can be deduced from observations of its gravitational 
radiation, at least if only the Newtonian orbit is important. Moreover, if 
one can measure the field amplitude (e.g. h'^'^^'^) plus fig^, and f2g«,, one 
can deduce from these the value of M and the distance r to the system! 
A chirping binary with a circular orbit, observed in gravitational waves, is 
a standard candle: one can infer its distance purely from the gravitational 
wave observations. To do this one needs the full amplitude, not just its 
projection on a single detector, so one generally needs a network of detec- 
tors or a long-duration observation with a single detector to get enough 
information. 

It is very unusual in astronomy to have standard candles, and they 
are highly prized. For example, one can in principle use this information 
to measure Hubble's constant. [26] 

6.1.1 Corrections 

In the calculation above we made several simplifying assumptions. For 
example, how good is the assumption that the orbit is circular? The Hulse- 
Taylor binary is in a highly eccentric orbit, and this turns out to enhance 
its gravitational wave luminosity by more than a factor of 10, since the 
elliptical orbit brings the two stars much nearer to one another for a period 
of time than a circular orbit with the same period would do. So there are 
big corrections for this system. 

However, systems emitting at frequencies observable from ground- 
based interferometers are probably well-approximated by circular orbits, 
because they have arrived at their very close separation by gravitational- 
wave-driven in-spiral. This process removes eccentricity from the orbit 
faster than it shrinks the orbital radius, so by they time they are observed 
they have insignificant eccentricity. 

Another assumption is that the orbit is well described by Newtonian 
theory. This is not a good assumption in most cases. Post-Newtonian orbit 
corrections will be very important in observations. This is not because the 
stars eventually approach each other closely. It is because they spend a 
long time at wide separations where the small post-Newtonian corrections 
accumulate systematically, eventually changing the phase of the orbit by an 
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observable amount. So it is very important for observations that we match 
signals with a template containing high-order post-Newtonian corrections, 
as described in Blanchet's lecture. But even so, the information contained 
in the Newtonian part of the radiation is still there, so all our conclusions 
above remain important. 



6.2 The r-modes 

We consider rotating stars in Newtonian gravity and look at the effect that 
the emission of gravitational radiation has on their oscillations. One might 
expect that the loss of energy to gravitational waves would damp out any 
perturbations, and indeed this is normally the case. However, it was a 
remarkable discovery of Chandrasekhar[3J that the opposite sometimes 
happens. 

A rotating star is idealized as an axially symmetric perfect-fluid sys- 
tem. In the Newtonian theory the pulsations of a perturbed fluid can be 
described by normal modes which are the solutions of perturbed Euler and 
gravitational field equations. If the star is stable, the eigenfrequencies a of 
the normal modes are real; if the star is unstable, there is at least one pair 
of complex-conjugate frequencies, one of which represents an exponentially 
growing mode and the other a decaying mode. (We take the convention 
that the time-dependence of a mode is exp{iat).) 

In general relativity, the situation is in principle the same, except 
that there is a boundary condition on the perturbation equations that 
insists that gravitational waves far away be outgoing, i.e. that the star loses 
energy to gravitational waves. This condition forces all eigenfrequencies to 
be complex. The sign of the imaginary part of the frequency determines 
stability or instability. 

The loss of energy to gravitional radiation can destabilize a star that 
would otherwise (i.e. in Newtonian theory) be stable. This is because it 
opens a pathway to lower-energy configurations that might not be acces- 
sible to the Newtonian star. This normally happens because gravitational 
radiation also carries away angular momentum, a quantity that is conserved 
in the Newtonian evolution of a perturbation. 

The sign of the angular momentum lost by the star is a critical di- 
agnostic for the instability. A wave that moves in the positive angular 
direction around a star will radiate positive angular momentum to infinity. 
A wave that moves in the opposite direction, as seen by an observer at 
rest far away, will radiate negative angular momentum. In a spherical star, 
both actions result in the damping of the perturbation because, for exam- 
ple, the postivive-going wave has intrinsically positive angular momentum, 
so when it radiates its angular momentum decreases and so its amplitude 
decreases. Similarly, the negative- going wave has negative angular mo- 
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mentum, so when it radiates negative angular momentum its amplitude 
decreases. 

The situation can be different in a rotating star, as first pointed out 
by Friedman and Schutz425| The angular momentum carried by a wave de- 
pends on its pattern angular velocity relative to the star 's angular velocity, 
not relative to an observer far away. If a wave pattern travels backwards 
relative to the star, it represents a small effective slowing down of the star 
and therefore carries negative angular momentum. This can lead to an 
anomalous situation: if a wave travels backwards relative to the star, but 
forwards relative to an inertial observer (because its angular velocity rel- 
ative to the star is smaller than the star's angular velocity), then it will 
have negative angular momentum but it will radiate positive angular mo- 
mentum. The result will be that its intrinsic angular momentum will get 
more negative, and its amplitude will grow. 

This is the mechanism of the Chandrasekar-Friedman-Schutz (CFS) in- 
stability. In an ideal star, it is always possible to find pressure-driven waves 
of short enough wavelength around the axis of symmetry (high enough an- 
gular eigenvalue m) that satisfy this condition. But it turns out that even 
a small amount of viscosity can damp out the instability in such waves , so 
it is not clear that pressure-driven waves will ever be significantly unstable 
in realistic stars. 

However, in 1997 Andersson[TJ pointed out that there was a class of 
modes called r-modes (Rossby modes) that no-one had previously investi- 
gated, and that were formally unstable in all rotating stars. Rossby waves 
are well-known in oceanography, where they play an important role in en- 
ergy transport around the Earth's oceans. They are hard to detect, having 
long wavelengths and very low density perturbations. They are mainly 
velocity perturbations of the oceans, whose restoring force is the Coriolis 
effect, and that is their character in neutron stars too. Because they have 
very small density perturbation, the gravitational radiation they emit is 
dominated by the current-quadrupole radiation. 

For a slowly-rotating, nearly-spherical Newtonian star, the following 
velocity perturbation is characteristic of r-modes: 

<5«" = ^(r)e'^''^V6rVcri„, (6.9) 

where <^{r) is some function of r determined by the mode equations. This 
velocity is a curl, so it is divergence- free; since it has no radial component, 
it does not change the density. If the star is perfectly spherical, these per- 
turbations are simply a small rotation of some of the fluid, and it continues 
to rotate. They have no oscillation, and have zero-frequency. 

If we consider a star with a small rotational angular velocity 51, then 
the frequency a is no longer exactly zero and a Newtonian calculation to 
first order in D, shows that there is a mode with pattern speed ujp = —a/ni 
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equal to 



uj„ — Q 



1 



1(1 + 1) 



(6.10) 



These modes are now oscillating currents that move (approximately) along 
the equipotential surfaces of the rotating star. 

For I ^ 2, ijjp is positive but slower than the speed of the star, so by the 
CFS mechanism these modes are unstable to the emission of gravitational 
radiation for an arbitrarily slowly rotating star. 

The velocity pattern given in Equation (16. 9p for {I — 2, m — 2) is 
closely related to the wheel model we described for current-quadrupole 
radiation in Figure 15.11 Take two such wheels and orient their axels along 
the X- and y-axes, with the star rotating about the z-axis. Choose the sense 
of rotation so that the wheels at positive- a; and positive-y are spinning in 
the opposite sense at any time, i.e. so that their adjacent edges are always 
moving in the same direction. Then this relationship will be reproduced 
for all other adjacent pairs of wheels: adjacent edges move together. 

When seen from above the equatorial plane, the line-of-sight momenta 
of the wheels reinforce each other, and we get the same kind of pattern that 
we saw when looking at one wheel from the side. However, in this case the 
pattern rotates with the angular velocity 2r2/3 of Equation (|6.10p . Since the 
pattern of line-of-sight momenta repeats itself every half rotation period, 
the gravitational waves are circularly polarized with frequency 4J7/3. Seen 
along the x-axis, the wheel along the x-axis contributes nothing, but the 
other wheel contributes fully, so the radiation amplitude in this direction 
is half that going out the rotation axis. Seen along a line at 45° to the 
a;-axis, the line-of-sight momenta of the wheels on the front part of the star 
cancel those at the back, so there is no radiation. Thus, along the equator 
there is a characteristic series of maxima and zeros, leading to a standard 
m = 2 radiation pattern. This pattern also rotates around the star, but 
the radiation in the equator remains linearly polarized because there is only 
the (g) component, not the ©. Again, the radiation frequency is twice the 
pattern speed because the radiation goes through a complete cycle in half 
a wave rotation period. 

This discussion cannot go into the depth required to understand the 
r-modes fully. There are many issues of principle: what happens beyond 
linear order in i7; what happens if the star is described in relativity and 
not Newtonian gravity; what is the relation between r-modes and the so- 
called 5-modes that can have similar frequencies; what happens when the 
amplitude grows large enough that the evolution is non-linear; what is the 
effect of magnetic fields on the evolution of the instability? The literature 
on r-modes is developing rapidly. We have included references where some 
of the most basic issues are discussed, [T71 [T51 [2S1 [201 ISI] but the interested 
student should consult the current literature carefully. 
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6.2.1 Linear growth of the r-modes 

We have seen how the r-mode becomes unstable when coupled to gravita- 
tional radiation, and now we turn to the practical question: is it important. 
This will depend on the balance between the growth rate of the mode due 
to relativistic effects and the damping due to viscosity. 

When coupled to gravitational radiation and viscosity, the mode has 
a complex frequency. If we define 5((t) := 1/r, then r is the characteristic 
damping time. When radiation and viscosity are treated as small effects, 
their contributions to the eigenfrequencies add, so we have that the total 
damping is given by 

111111 , , 

^7^ = + -, - = - + -, (6.11) 

r(i2) TGR Tv Tv Ts Tb 

where I/tgr, 1/tv are the contributions due to gravitational radiation emis- 
sion and viscosity, and where the latter has been further divided between 
shear viscosity (1/ts) and bulk viscosity (1/Tb). 

If we consider a "typicaf neutron star with a polytropic equation of 
state p = kp^ (for which k has been chosen so that a 1.5 Mq model has a 
radius R — 12.47 km), and if we express the angular velocity in terms of 
the scale for the approximate maximum speed \pKGp and the temperature 
in terms of 10^ K, then it can be shown thatj30j 



1 . _L ('i?!iy% _L rii^^rviav + -L r^.V , (6,2) 



P j T^\P I \ T j T„ \\tS>Y.j 



where the scaling parameters fgv, fbv, Tgw and the exponents Pgw and pbv 
have to be calculated numerically. Some representative values relevant to 
the r-modes with 2 ^ ^ ^ 6 are in Table lO [30] . 

Table 6.1. Gravitational radiation and viscous time scales, in seconds. Negative 
values indicate instability, i.e. a growing rather than damping mode. 



I 


m 


T-9™(S) 


Pgw 


Tbv (s) 


Pbv 


Tsv (s) 


2 
3 


2 
3 


-20.83 
-316.1 


5.93 
7.98 


9.3 X 10^° 
1.89 X 10^" 


1.77 
1.83 


2.25 X 10** 
3.53 X 10^ 



The physics of the viscosity is interesting. It is clear from Equation 
(J6.12I) that gravitational radiation becomes a stronger and stronger desta- 
bilizing influence as the angular velocity of a star increases, but the vis- 
cosity is much more complicated. There are two contributions: shear and 
bulk. Shear viscosity comes mainly from electrons scattering off protons 
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and other electrons. This effect fahs with increasing temperature, just as 
does viscosity of every-day materials. So a cold, slowly rotating star will 
not have the instability, where a hotter star might. But at high tempera- 
tures, bulk viscosity becomes dominant. This effect arises in neutron stars 
from the nuclear physics. Neutron-star matter always contains some pro- 
tons and electrons. When it is compressed, some of these react to form 
neutrons, emitting a neutrino. When it is expanded, some of the neutrons 
beta-decay to protons and electrons, again emitting a neutrino. The emit- 
ted neutrino is not trapped in the star; within a short time, of the order of 
a second or less, it escapes. This irreversible loss of energy each time the 
star is compressed creates a bulk viscosity. Now, bulk viscosity acts only 
due to the density perturbation, which is small in r-modes. So the effect 
of bulk viscosity only dominates at very high temperatures. 

The balance of the viscous and gravitational effects is illustrated in 
Figure [B.1J 30'. This is indicative, but not definitive: much more work is 
needed on the physics of viscosity and the structure of the modes at large 
values of ft (small P) . 

6.2.2 Nonlinear evolution of the star 

Our description so far is only a linear approximation. To understand the full 
evolution of the r-modes we have to treat the non-linear hydrodynamical 
effects that become important as the modes grow. This could only be done 
with a numerical simulation, which some groups are now working on. But 
it is possible to make simple estimates analytically. 

We characterize the initial configuration with just two paramters: the 
uniform angular velocity J7, and the amplitude a of the r-modes perturba- 
tion. The star is assumed to cool at the accepted cooling rate for neutron 
stars, independently of whether it is affected by the r-mode instability or 
not. The star is assumed to lose angular momentum to gravitational radia- 
tion at a rate given by the linear radiation field, with its large amplitude a. 
This loss is taken to drive the star through a sequence of equilibrium states 
of lower and lower angular momentum. Details of this approximation are 
in|31). here we report only the results. The evolution turns out to have 
three phases: 

i) Initially the angular velocity il of the hot rapidly rotating neutron 
stars is nearly constant, evolving on the viscous time-scale I/tv, while 
the amplitude a grows esponentially on the gravitational radiation 
time-scale I/tgr. 

ii) After a short time non-linear effects become important and stop the 
growth of the amplitude a. Most of the initial angular momentum of 
the star is radiated away by gravitational radiation. The star spins 
down and evolves to a point where the angular velocity $7 and the 
temperature is sufRciently low that the r-mode is stable. 
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Figure 6.1. The balance of viscous and gravitational radiation effects in the 
r-modcs is illustrated in a diagram of rotation speed, showing the ratio of the 
maximum period Pk to the rotation period P, versus the temperature of the 
star. The solid curve indicates the boundary between viscosity-dominated and 
radiation-dominated behavior: stars above the line are unstable. The dashed 
curves illustrate possible nonlinear evolution histories as a young neutron star 
cools. 




10 11 



iog,oT/iK 



Conclusion 69 

iii) Finally gravitational radiation and viscosity damp out the r-niode and 
drive the star into its final equilibrium configuration. 

This may take about a year, a timescale governed by the cooling time 
of the star. During this year, the star would radiate away most of its angu- 
lar momentum and rotational kinetic energy. This could be a substantial 
fraction of a solar mass in energy. 

6.2.3 Detection of r-mode radiation 

The large amount of energy radiated into the r-modes makes them at- 
tractive for detection, but detection will not be trivial. The r-mode event 
occurs at the rate of supernovae: some fraction (hopefully large) of all 
supernovae leave behind a rapidly spinning neutron star that spins down 
over a 1-year period. This means we should have sufficient sensitivity to 
reach to the Virgo Cluster (20 Mpc distance). Estimates [21] suggest that 
a neutron star in the Virgo Cluster could be detected by second genera- 
tion of LIGO and VIRGO gravitational wave detectors with an amplitude 
signal-to- noise of about 8, provided one can use matched filtering (exact 
template matching). 

It will not be easy to use matched filtering, since one must follow all 
cycles of the signal as the star spins down, and we won't know this well 
because of many uncertainties: initial temperature, initial spin distribution, 
detailed physics of viscosity, and so on. But it would be helpful to have a 
parametrized model to take account of the uncertainties, so that we could 
look for a significant fit to one or more of the parameters. 

In addition, it is likely that, if a significant proportion of all neutron 
stars went through the r-mode instability, then the Universe has been filled 
by their radiation. There should be a background with an energy density 
flgu) that is a good fraction of the closure density. Its lower frequency 
limit should be around 200 Hz in the rest frame of the star. When we see 
radiation cosmologically, its lower frquency limit will indicate the epoch at 
which star formation began. 

It is clear that the discovery of this new source of gravitational waves 
will open several prospects for astronomy. Observations could be used as 
supernovae detectors, revealing supernovae hidden in clouds of dust, iden- 
tifying them about at year after they are formed. The existence of the 
radiation raises several prospects and questions about the physics of neu- 
tron stars, not least the interaction of magnetic fields with the instability. 

6.3 Conclusion 

These lectures have taken us through the basic theory of gravitational ra- 
diation and its applications in astrophysics, so far as we can understand 
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and predict them now. In a few years, perhaps as httle as 2, perhaps as 
many as 8, we wih start to make observations of gravitational radiation 
from astrophysical sources. If gravitational wave astronomy follows other 
branches of observational astronomy, it will not be long before completely 
unexpected signals are seen, or unexpected features in long-predicted sig- 
nals. To interpret these will require a joining a physical understanding of 
the relationship between gravitational radiation and its source to a wide 
knowledge of astronomical phenomena. I encourage the students who have 
attended these lectures, and others who may study them, to get themselves 
ready to contribute to this activity. It will be an exciting time! 
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Lecture 1 Exercises 

Exercise 1 

(a) Let us take the form of the wave to be 

h •' — e'^ h+{t — n ■ ii.) 

where e^ is the polarization tensor for the polarization, and where n is the 
unit vector in the direction of travel of the wave. We will let h+ be an arbitrary 
function of its phase argument. 

If the wave travels in the x-z plane at an angle 9 to the z-direction, then 
the unit vector in our coordinates is 

n* = (sin^, 0, cos 6). 

We need to calculate the polarization tensor's components in the x, y, z coor- 
dinates. We do this by rotating the © polarization tensor from its TT-form in 
coordinates parallel to the wavefront to its form in our coordinates. This requires 
a simple rotation around the j/-axis. The transformation matrix is: 

/ cos9 sine* \ 
A^'fe = 1 . 

\ -sine* cos6' / 

The polariation tensor in our coordinates (primed indices) becomes: 

e-' = A-' (A me 

(cos^ 6 — sin 6* cos 61 \ 

0-10 
-sin 6* cos 6* sin'' 9 ) 

Notice that the new polarization tensor is again traceless. 

The gravitational wave will be, at an arbitrary time t and position (a;, z) in 
our {x, 2:)-plane, 

h •' — e' h-t-{t — xsm9 — zcosff). 
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For this problem we need the a;2;-component because the photon is propagating 
along this direction, and we will always stay at 2: = 0, so we have 

h ^^ = cos 9h+(t — X sin 0). 

We see that for this geometry the wave amplitude is reduced by a factor of cos^ 9. 
Generalizing the argument in the text, the relation between time and posi- 
tion for the photon on its trip outwards along the i-direction ist = to + x, where 
to is the starting time. The analogous relation after the photon is reflected is 
t = to + L + {L — x), since in this case x decreases in time from L to 0. If we put 
these into the equation for the linearized corrections to the return time, we get 

treturn = to + 2L + - cos^ 9 I h+[to + (1 - sin 6')x]d2; 



/ /i+[to + 2L- (l + sin6')j;]da; 

^0 



This expression must be differentiated with respect to to to find the variation 
of the return time as a function of the start time. The key point is how to 
handle differentiation within the integrals. Consider, for example, the function 
h+[to + (1 — sin9)x]. It is a function of a single argument, 

C ~ to + (1 -sin6')i 

so derivatives with respect to to can be converted to derivatives with respect to 
X as follows 

dh+ dh+ d^ dh+ 

dto d^ dto df 

dh+ _ dh+ d^ _ , . a\'^^+ 

dx d^ dx d^ 

It follows that 

§i = fi/(l-sin.). 
dto dx 

On the return trip the factor will be —(1 + sinS). So when we differentiate we 

can convert the derivatives with respect to to inside the integrals into derivatives 

with respect to x. Taking account of the factor cos^ ^ = (I — sin^)(l + sin 6*) in 

front of the integrals, the result is 

%^ = 1 + i(l + sine) /^ ^[to + (1 - sm9)x]dx + 
dto 2 Jo dx 

l{l~sm9) [ ^[to + 2L-{1 + sin 9)x]dx. 



2 Jq dx 



The integrals can now be done, since they simply invert the differentiation by 
X. Evaluating the integrands at the end points of the integrals gives Equation 



'^^"'""" = 1 - -{l + sin9)h+{to) + sin9h+[to + (1 - sin6»)i] 
dto 2 

+ -(f-sin6')/i+(to + 2L). 
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(b) If we Taylor-expand this equation in powers of L about L = 0, the 
leading term vanishes, and the first-order term is: 

^hpIIL =Lsme{l - sme)h+{to) + 
dto 

L(l -sin6l)/i+(io), 

= Lcos 9h+{to). 

This is just what was required. The factor of cos^ 6 comes, as we saw above, from 
the projection of the TT field on the a;-coordinate direction. 

(c) All the terms cancel and there is no effect on the return time. 

Exercise 2 

This is part of the calculation in the previous example. All we need is the 
segment where the light travels from the distant end to the center: 

1 2 f^ 

tout = to -f - COS 9 h+[to + L — {1 + sin d)x]dx 

^ Jo 



and so %^ is: 

dtn 



%i = 1 + i(l - sine) [h+{to - sm9L) - h+{to + L)] 
dto ^ 



Exercise 3 

This question is frequently asked, but not by people who have done the 
calculation. The answer is that the two effects occur in different gauges, not in 
the same one. So they cannot cancel. The apparent speed of light changes in the 
TT-gauge, but then the positions of the ends remain fixed, so that the effect is 
all in the coordinate speed. In a local Lorentz frame tied to one mass, the ends 
do move back and forth, but then the speed of light is invariant. 

Exercise 4 

To the first order we have 

^afiv ^ '^ av,fi '^cifi,v^ 
^ al3,i^ ~ nV {hap^av + t>'CTa,Pv — hal3,<yv) , (i) 

'^av,fi ~ 9*^ {hau,aP + haa,ufi — hai',al3) ■ (h) 

The gauge transformation for a perturbation in linearized theory is 

haf) = hafi — ia,[i — £,l3,a- (hi) 

Substituting (iii) into (i) and (ii), we obtain 
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If we find the difference between the two formulas above we get 

Lecture 4 Exercises 

Exercise 5 

The action principle is: 



SI = hnvd X = — I G^" \/~^hnvd X = Q. (i) 

J t>5Mi' J 

If we perform a infinitesimal coordinate transformation x^ — > x^ + ^^ without 
otherwise varying the metric, then the action I must not change: 

2 j G^-'i^.^d^. 
This can be transformed in the following way: 

The first integral is a divergence and vanishes. The second, because of the arbi- 
trariness of ^;j, gives the Bianchi's identities: 

G^"., = 0. 



Exercise 6 

The two polarization components are /i+^ = — ft^" = A+e"**^*"^' and /i^" 
^^g-jfe(t-z) rpj^g energy flux is the negative of 



